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1. Introduction 

The nonlinear steepest descent method is an approach for determining the asymptotic behav¬ 
ior of solutions of matrix Riemann-Hilbert (RH) problems that depend on a large parameter. 
By deforming the contour in such a way that the jump is small everywhere except near a finite 
number of critical points, detailed asymptotic formulas can be obtained by summing up the 
contributions from the individual critical points. As the inverse scattering transform formalism 
expresses the solution of a nonlinear integrable PDE in terms of the solution of a RH problem, 
the nonlinear steepest descent method is particularly well-suited for determining asymptotic 
formulas for integrable equations. In fact, the method was first developed by Deift and Zhou in 
[TO] . where the long time behavior of the solution of the modified Korteweg-de Vries (mKdV) 
equation on the line was established. Given initial data u o(cc), x € R, in the Schwartz class, they 
showed that the asymptotics of the associated solution u(x, t ) can be determined by performing 
a triangular factorization of the jump matrix followed by a contour deformation. The latter 
step could be completed thanks to the introduction of an ingenious analytic approximation of 
the jump matrix. 


E-mail address : jlenells@kth.se. 
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In this paper, we consider the derivation of asymptotic formulas for initial- boundary value 
(IBV) problems. Our goal is to provide a rigorous and yet accessible treatment. We prove two 
nonlinear steepest descent theorems suitable for finding the asymptotics in the similarity and 
the self-similar sectors respectively. We then apply these theorems to derive uniformly valid 
asymptotic formulas. For definiteness, we consider the mKdV equation 

Ut “t“ Q'U U x U X xx — 0; (1-1) 

posed in the quarter plane {x > 0, t > 0}. 

We first consider the asymptotic behavior in the similarity sector defined by 

t > 1, x 3 /t ->oo, x < Nt , N constant. (1.2) 

The condition x 3 /t —> oo defines the left boundary of the sector in the (cc,t)-plane while the 
inequality x < Nt defines its right boundary. The solution u(x , t ) can be expressed [3] in terms 
of the solution M(x,t,k) of a RH problem whose jump matrix involves two spectral functions 
r(fc) and h(k) which are defined in terms of the initial data uo(x ) = u(x, 0), x > 0, and of the 
boundary values 

g 0 (t) = u(0,t), 5 i(t) = u x (0,t), g 2 (t) = u xx (0,t), t> 0. (1.3) 


In addition to a jump across the real axis, which is present also in the case of the initial 
value problem, the RH problem for the IBV problem also has jumps across the two lines Me 5 ? 
and Ke^, see Figure [lj The jumps across these two lines involve the spectral function h(k), 
whereas the jump across R involves the spectral function r(k). There are two critical points 


located at ±ko where fco = yj yfj- At a formal level, it is relatively straightforward to obtain 
the leading order asymptotics of u(x , t) by reading off the contributions from the two critical 
points. However, the rigorous derivation of a uniform expansion with precise error terms is more 
involved. We establish such an expansion (see Theorem 4.1l by implementing the following seven 
steps: 


1. Find an analytic approximation h a (t,k ) of h(k). 

2. Deform the contour to eliminate the part of the jump that involves h a (t , k). 

3. Conjugate the RH problem to arrive at an appropriate triangular factorization. 

4. Find an analytic approximation of the jump matrix along R. 

5. Deform the contour so that the jump is small everywhere except near ±fco- 

6. Apply a nonlinear steepest descent theorem to find the asymptotics of M(x,t , k). 

7. Find the asymptotics of u(x,t). 


Steps 3-7 have analogs in the case of the initial value problem, whereas steps 1 and 2 are unique 
to the half-line problem. 

We also consider the asymptotic behavior of u{x, t) in the self-similar sector S defined by 


S = {(a;,f) 1 1 > 1 andO < x < Nt 1// ' 3 }, N constant. (1.4) 


In this case, the two critical points ±fco merge at the origin as t —> oo. Through a series of 
seven steps, conceptually similar to the seven steps listed above, we derive (see Theorem 7.1) 
an asymptotic formula for the solution. 

Our derivations require that r(k) and h{k ) possess a certain amount of regularity. Since 
regularity of the spectral functions corresponds to decay of the initial and boundary values, 
this means that the asymptotic formulas only are valid provided that the boundary values in 
(1.3) have some decay as t —> oo (see Theorem 2.1 and Remark 4.2 for details). We will not 
consider the interesting, but difficult, question of deriving asymptotic formulas for problems 
whose boundary values do not vanish for large t. 

Other works studying asymptotics in the context of IBV problems for nonlinear integrable 
PDEs include GO E ED IH IB]. It was shown in [3] that the solution of the mKdV equation ( |1.1[ ) 
on the half-line is 0(f -1,/2 ) in the asymptotic sector 0 < c\ < x/t < C 2 . Our asymptotic formula 
in the similarity region is stronger than this result, because it provides an exact expression for 
the coefficient of f -1 / 2 in the asymptotic expansion. In mm the solution of the defocusing 
nonlinear Schrodinger (NLS) equation on the half-line was shown to be 0(t -1 / 2 ) for 0 < c\ < 
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x/t < C 2 , and, for the focusing NLS, formulas for the asymptotic solitons were presented. Similar 
results were obtained for the KdV equation on the half-line in m The asymptotic behavior of 
the solution of the Camassa-Holm equation on the half-line was analyzed in [8]. 

1.1. Main results. The paper presents four main theorems of two different types: 


(z) Theorem 3.1 and Theorem 6.1 are nonlinear steepest descent theorems. They are of a 


more general character and are suitable for determining asymptotics of two different 
classes of RH problems depending on a large parameter. Although they are formulated 


with the quarter plane problem for the defocusing mKdV equation (1.1) in mind, they 


could also be relevant for other initial-boundary value problems (perhaps after appropri¬ 
ate adjustments for different symmetries or dispersion relations). Indeed, the goal has 
been to formulate these theorems assuming only those properties of the RH problem 


associated with (1.1) which are relevant for the application of the nonlinear steepest 
descent method. 


(zz) Theorem 4.1 and Theorem 7.1 are asymptotic theorems specific to the quarter plane 


problem for equation (1.1). Theorem 4.1 establishes the asymptotics of the solution in 


the similarity sector and is an application of Theorem |3.1| t o the RH problem associated 
with <o>- Similarly, Theorem |7.1| (see also Corollary |7.2[ ) establishes the asymptotics 
of the solution in the self-similar sector and is an application of Theorem |6.1| to the same 
RH problem. 

We hope that this division of the derivation into two steps will make it easier when applying 
the approach to other initial-boundary value problems in the future. We also hope that it 
clarifies the logical structure of the proofs by isolating those properties of the RH problem that 
are essential for the nonlinear steepest descent arguments. 

1.2. Effect of the boundary. The asymptotic formulas for the half-line problem obtained in 
Theorems |4.1| and |7.1| are similar to the analogous formulas for the pure initial value problem. 
In fact, in the similarity sector, the asymptotic formulas on the line and on the half-line have the 
exact same functional form. The only difference is that for the half-line problem, the definition 
of the spectral function r{k), which enters the asymptotic formulas, in addition to the initial 
data, also involves the boundary values (see Remark 4.3). In other words, the only effect of 


the boundary is to modify r(k). Intuitively, we can understand this as follows: Since x grows 
faster than t 1 ' 3 in the similarity region, the distance to the boundary eventually gets so big that 
what happens at the boundary has a small effect on the solution (recall that we assume that 
the boundary values decay as t — > oo). This means that, from the point of view of the leading 
order asymptotics in the simlarity sector, the boundary values and the initial data play similar 
roles. It is therefore not surprising that they enter the asymptotic formula in similar ways. 

In the self-similar sector S, the effect of the boundary is much greater. For the problem on 
the line, the asymptotics of the solution u(x, t) in this region is given by 


fix, t) = 


’( 


( 3 t )!/ 3 : 


5 , 0 ,- s ) 


(3 t)i 


+ 0(t 3 ), t —> oo, (x,t) £ 


(1.5) 


where 


u* denotes a particular solution of the Painleve II equation specified by a complex 
parameter s £ zR, |s| < 1 (see Appendix [A| . For the half-line problem, the asymptotic analysis 
of the relevant RH problem leads to the same asymptotic formula, except that the value of the 
parameter s now also depends on the boundary values in addition to the initial data. So far, 
the situation is analogous to what happens in the similarity region. However, if we consider a 


quarter plane solution of (1.1) with sufficient smoothness and spatial decay, then the associated 


spectral functions are not independent, but are related by the so-called global relation. It turns 
out that, since we are assuming that the boundary values decay as t —> oo, the global relation 
forces the parameter s to vanish. The solution u p corresponding to s = 0 vanishes identically, 
so this means that, for the half-line problem, we have u(x,t ) = 0{t~ S) as t —> oo in the self- 
similar region. This observation can be interpreted as follows: In the self-similar sector, since 
x < Nt 1 / 3 , the boundary is close enough to have a noticeable effect on the solution. Since we 
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assume that the boundary values decay to zero as t —> oo, the solution gets drawn down to zero 
faster in the presence of the boundary than in its absence. 


1.3. Organisation of the paper. In Section[2j we review the construction of solutions of (1.1) 
in the quarter plane via RH techniques. In Sections [3] and [6j we prove the nonlinear steepest 
descent results Theorem |3.1| and Theorem |6.1[ respectively. In Sectio ns [4] and [7| we apply these 
theorems to find the asymptotics of the quarter plane solution of (1.1) in the similarity and 
self-similar sectors, respectively. A few facts related to the RH problem associated with the 
Painleve II equation are collected in the appendix. These facts are used in Section [5] to establish 
the solution of a model RH problem which is essential for the proof of Theorem |6.1| 


1.4. Open problems. We end the introduction by listing three interesting, largely open, prob¬ 
lems: 


1. In this paper we determine the asymptotic behavior of u(x 7 1) provided that all the boundary 
values {ft£zi(0, t)}^ are known. However, for a well-posed problem, only a subset of the initial 
and boundary values can be independently prescribed. If all boundary values are not known, 
our asymptotic formulas (see Theorem 4.1 and Theorem 7.1) still provide some information 
on the solution, but since the function r(k) is unknown, the exact form of the asymptotics 
remains undetermined. It is an outstanding open problem to derive asymptotic formulas 
which involve only the prescribed data. 

2. The asymptotic formulas of Theorem |4.1| and Theorem |7.1| only apply provided that the 
boundary values gj{t) = d^.u(0,t), j = 0,1, 2, possess a certain amount of decay as t —> oo. 
For example, in view of Theorem |2.1| the following decay assumption is sufficient: 

(l + t) 11 fl? ) (*)e£ 1 ([0,c»)), j = 0,1,2, * = 0,1,... ,4. 


However, as already mentioned, only a subset of the boundary values can be prescribed for 
a well-posed problem. Thus, even though the given boundary data decay as t —> oo, it is 
not immediate that all the boundary values decay as t —> oo. Is it true that decaying data 
leads to decaying boundary values in general? The investigation of this question involves an 
analysis of the generalized Dirichlet to Neumann map. Important progress in this direction 
has been made by Antonopoulou and Kamvissis for the NLS equation on the half-line. For 
the defocusing NLS equation with vanishing initial data, they have shown that if the Dirichlet 
data decay as t —> oo, then so does the Neumann value pQ. 

3. Another challenging but very important problem consists of deriving asymptotic formulas 
for solutions whose boundary values do not decay as t —► oo. In this case, the definition 
of the spectral functions r(k) and h{k) has to be modified. If the asymptotic behavior of 
the boundary values as t —> oo is known, substitutes for r(k) and h{k) can be defined by 
subtracting off the known asymptotics in the relevant integral equation. Of particular interest 
in this regard are problems with (asymptotically) time-periodic data. For the NLS equation 
with periodic boundary conditions, some pioneering results in this direction can be found in 

m\- 


2. Quarter plane solutions 


The rnKdV equation (1.1) admits the Lax pair 

ik[a 3 ,n] = U/t, 

|+4ifc 3 [(7 3 ,/t] = V/t, 

where g{x, t, k ) is a 2 x 2-matrix valued eigenfunction, k £ C is the spectral parameter, and 


( 2 . 1 ) 


C 3 


! 0 

0 -1 


li(x,t) = 


0 u 
u 0 


V(x, t , k) = 


—2iu 2 k 

—Auk 2 — 2 iu x k — 2 u 3 + u xx 


—Auk 2 + 2 iu x k — 2 u 3 + u x 
2 iu 2 k 
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Figure 1. The contour E and the domains {Djin the complex k-plane. 


We define the spectral functions {a(k),b(k),A(k) 1 B(k)} by 


X(0,fc) 


(aW K k )\ 

\b(k) a(k)J ’ 


T(0, k) 


A{k) B(k)\ 

B(k) A{k) J ’ 


( 2 . 2 ) 


where X(x, k) and T (t, k) are the solutions of the Volterra integral equations 


X(x, k) = I + 
T(t,k) =1 + 


a -ik(x'-x)a 3 | | ( I 


U(a/, 0)X(a/, k)dx ’, 

e 4ik 3 {t'-t)a 3 V ( 0) t / ) jfe)^ 


and e CT3 acts on a 2 x 2 matrix M by e a3 M = e" 3 Me °’ 3 . The open domains {Dj}{ of the 
complex fc-plane are defined by (see Figure [l]) 

Di = {Imfc < 0} n {Imfc 3 > 0}, D 2 = {Im/c < 0} n {Im/c 3 < 0}, 

D 3 = {Imfc > 0} n {Imfc 3 > 0}, Di = {Im/c > 0} n {Imfc 3 < 0}. 


We let E denote the contour separating the Dj 's oriented as in Figure |Tj 
We define h(k) and r(k) by 


w = -jmm- k e (2 - 3 “» 

r {k) = ^\ + h(k), k G R, (2.3b) 

a{k) 

where d{k) = a(k)A(k) — b(k)B(k) and D 2 denotes the closure in C of the set D 2 . It was 
shown in HZ! that if the initial data uq(x) = u(x, 0) and the boundary values gjit) = d{,u(0,t), 
j = 0,1,2, satisfy 


r« 0 GC TO + 1 ([0,oo)), 

I (1 -\-x) n u^\x) e T^fOjOo)), i = 0, 

^^([0,00)), j = 0,1,2, [ ~ A) 

l(H-t) n 5f(t)Gi 1 ([0,oo)) ) i = 0,1,2, * = 0,1,..., [=*£=<], 

for some integers n,m > 1, and the spectral functions satisfy the so-called global relation 

A(k)b(k) — B{k)a(k) = 0, k £ D 3 , (2-5) 


then r(k) and h{k ) have the following properties: 
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• re C n (R). 

• h(k) is analytic in D 2 and h^\k) has a continuous extension to D 2 for each j = 
0,1,... ,n. 

• There exist complex constants {hi}™ such that, for each j = 0,1,..., n, 

r«(fc) =0(r m - 1+2i ), |A:| — »• 00, fceR, 


s?(T + - + £) + 0(t 


—m— l+ 2 j 


k e d 2 . 


( 2 . 6 ) 

(2.7) 


• r(k) = r{—k) for k e R and h(k) = h(—k ) for k e i? 2 - 

• sup fe 6 R |r(fc)| < 1 . 

The following theorem expresses the quarter plane solution of 0 in terms of the solution 
of a 2 x 2-matrix RH problem (see Theorem 7 of JT7 for a detailed proof). 


Theorem 2.1. Suppose ug, go, ffi) satisfy (2-4) with n = 1 and m = 4. Define the jump 
matrix J(x,t,k) by 


J(x, t,k) = < 


1 ° 

h{k)e~ 2ikx+8ikH l)' 

\ _ r ^J.'j e 2ikx—8ik 3 * D 

r(k) e - 2ikx+&ikH 1 — |r(fc )| 2 y 
1 h(fi)e 2ikx - 8ik3t \ 

.0 1 J ’ 


k £ dDi, 
k £ R, 
k £ dD 4 . 


( 2 . 8 ) 


Suppose the homogeneous RH problem determined by ( Y,,J{x,t , •)) /ias only the trivial solution 
for each (x,t) £ [0, 00 ) x [0, 00 ). Suppose the spectral functions satisfy the global relation (2.5). 
Then the L 2 -RH pro61er7|^] 


j M(x,t,-) £l + E 2 ( C\E), 

|M + (a;, t, k) = M_(x, t, k)J(x, t, k) for a.e. k £ E, 
has a unique solution for each ( x,t ) £ [0, 00 ) x [0,oo). Moreover, the nontangential 


(2.9) 


u(x,t) = —2i lim (kM(x,t,k)) 12 


k—t 00 


( 2 . 10 ) 


exists for each (x,t) £ [0, 00 ) x [0, 00 ) and the function u(x,t) defined by (2.10) is a solution of 
the mKdV equation 0 in the quarter plane {x > 0 , t > 0 } such that 

(a) u : [ 0 , 00 ) x [ 0 , 00 ) —)• R is C 8 in x and C 1 in t. 

(b) u(0,t) = g 0 (t), u x (Q,t) = 3 i(t), and u xx (0,t) = g 2 (t) for t > 0 . 

(c) u(x, 0 ) = uo(x) for x > 0 . 

The long time asymptotics of the quarter plane solution u{x,t) can be determined via a 
nonlinear steepest descent analysis of the RH problem (2.9). The next section contains a non¬ 
linear steepest descent theorem (Theorem 3.1) suitable for finding the asymptotics of u in the 
similarity sector ( 1 . 2 ). 


Remark 2.2. Equation (2.6) shows that for smooth initial and boundary values, the spectral 

00 . Thus, as in the case of the initial value problem, the 
A 00 . In contrast, regardless of how smooth the initial and 


function r(k) decays quickly as k 
jump across R vanishes quickly as k 


boundary values are, the function h(k) in general only decays as l/k as k —> 00 , see (2.7). We 
will overcome this difficulty by making a careful choice of the analytic approximation h a . 


1 FI ere the generalized Smirnoff class E 2 (C \ E) consists of functions f(k) analytic in C \ E with the property 
that for each component Dj of C \ E there exist curves {C n }J° in Dj such that C n eventually surround each 
compact subset of Dj and sup n>1 | /j| j < 00. 

2 X 

The notation lim indicates the limit as k £ C approaches 00 nontangentially with respect to T. 

k—too 

















THE NONLINEAR STEEPEST DESCENT METHOD 


7 



Figure 2 . The contour X = X\ U • • • U X 4 . 


3. A NONLINEAR STEEPEST DESCENT THEOREM 


The goal of this section is to prove a nonlinear steepest descent result (Theorem |3.1| ) suitable 
for determining asymptotics of a class of RH problems which arise in the study of long-time 
asymptotics in the similarity sector. Although the theorem is formulated with equation in 
mind, analogous results can be formulated for other IBV problems. Before stating the theorem, 
we describe the set-up and make a number of assumptions (see Remark 3.2 for a short discussion 
of the intuition behind these assumptions). 

Let X denote the cross X = X\ U • • • U X 4 C C where the rays 


Xi = {seT | 0 < s < 00 }, X 2 = {seT | 0 < s < 00 }, 

X3 = {se __r | 0 < s < 00 }, X4 = {se _T | 0 < s < 00 }, (3.11) 

are oriented as in Figure [2] For r > 0, let X r = X[ U • • • U X 4 denote the restriction of X to 

the open disk of radius r centered at the origin, i.e. X r = X C {\z\ < r}. 

Let 1 C K be an interval. Let p,e,k 0 : I —> (0, 00 ) be bounded strictly positive functions 
such that e(£) < fco(£) f° r C £ We henceforth drop the £ dependence of these functions and 
write simply p, e, fco for p(£), e(£), &o(£)> respectively. 

Let r = T(£) be a family of oriented contours parametrized by £ € I and let 


f := r U {k £ C | \k ± k 0 \ = e} 


denote the union of T with the circles of radius e centered at ±fco oriented counterclockwise. 
Assume that, for each £ £ I: 

(n) F and T are Carleson jump contours up to reorientation of a subcontour. 

(T2) r contains the two crosses ±fco + as oriented subcontours. 

(T3) F is invariant as a set under the map k 1 — > —k. Moreover, the orientation of T is such that 
if k traverses T in the positive direction, then —k traverses T in the negative direction. 
(T4) The point 00 in the Riemann sphere can be approached nontangentially with respect to 

F. 


Assume that the Cauchy singular operator Sp defined by 


( Sph)(z ) = lim — [ ) d z ' 

r->0 7 Tl Jf\{z' | \z'-z\<r} Z ~~ Z 


(3.12) 


is uniformly bounded on L 2 (T), i.e. sup^ eI ll^f llB(L 2 (f)) < °°- 

Consider the following family of L 2 - RH problems parametrized by the two parameters £ £ X 
and t > 0: 


m(£,f, •) £ I + E 2 (C\T), 

m+ (£, t, k ) = m_ (£, f, k)v{C >1 1 , k) for a.e. k £ T, 


(3.13) 
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where the jump matrix v((, t, k) satisfies 

»((, t, •) := r(c,i, •) - / e i^r) n i 2 (r) n i°°(r), Cez, t > o, (3.14) 


det v((, t, •) = 1 a.e. on T, ( G Z, t > 0, 


and 


v((,t,k) = v((,t,—k), C € X, t > 0, k £ F. 


(3.15) 

(3.16) 


Let t := tp 2 . Let T.y be the union of the two crosses ±ko + X e and let T' = r\r.Y- Suppose 

r —> oo, (€l, (3.17) 

uniformly with respect to £ £ Z. Moreover, suppose that near kg the normalized jump matrix 


KC>V)IUi(r') = 0(er *), 
KC,V)IU~(r') = 0(0, 


vo{(,t,z) = (t 3 u C,t,h 0 -- ct 3 , z £ X p , 

P 


(3.18) 


has the form 


v 0 (C,t,z) = C 


1 0 ] 

R l {C,,t,z)z~ 2 iv ^e t ^ z '> l)' 

1 -~R 2 (C,t,z)z 2i ^e- t<l> ^’ z '> S 

.0 1 

1 O' 

^-R 3 (C,t,z)z- 2 iv ^'>e t ‘ l> ^ z '> 1 / 

1 R 4 {( 1 t,z)z 2i, '^e- t ^ c - 2 ) n 

.0 1 


z £ X p , 
Z £ X£, 


Z £ Xg, 


z £ X%, 


(3.19) 


where: 


The phase (/)((, z ) is a smooth function of ((,z)€lxC such that 


and 


4>((, 0 ) £ iK, 

J 

g(CO)-o, 

g(C.O)=., <€l, 

(3.20) 

\z \ 2 

Re<j)(C,z) < - — , 


zgx fuY|, <£i, 

(3.21a) 

\z\ 2 

Re (/>((, z) > —, 


z£X p ux p , c eX, 

(3.21b) 

*(C,*)-#C,o)-^ 

<cliC, 

p 

zei", C ex, 

(3.21c) 


where C > 0 is a constant. 

• There exist a function q : X —> C with sup^ eI |g(C)| < 1, and constants (a,L) £ 
[|,1) x (0, oo) such that the functions {Rj((,t, z)}f satisfy the following inequalities: 


Ri(C,t,z) 

-q(.Q\<m a e^, 

z £ X p , 

R2(C,t,z) 

9(0 

T \ z \a tpr 

< X - e 6 

— 1 p \ 

z £ X %, 

l-|9(C)h 

Rs(C,t,z) 

9(C) 

i ict *M 2 

< L - e 6 

— I p \ 

z € -^ 3 > 

l-|9(C)r 

\Ri(C,t,z) 

-9(01 <Lf a e tJ ^, 

-2 € X^, 


() £ X, t > 0. 


(3.22) 


The function z/(£) is defined by v{Q = — X- ln(l — |g(C)| 2 )- 
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Theorem 3.1 (Nonlinear steepest descent). Under the above assumptions, the L 2 -RH pi'oblem 


(3.13) has a unique solution for all sufficiently large t and this solution satisfies 

i+a\ 


lim (km((,t,k )) 12 = _ 2 ‘ leRc ^(,t) + o( er - 


oo, ( 


k—> OO " ’ ' yj T 

where the error term is uniform with respect and (3{C,,t) is defined by 

f}(C,t) = )e®(« _ars 9 ^ +arsr ^ 1 '^^)e _t ^^ ,0 ^ _<1/ ^. 


(3.23) 


(3.24) 


The existence of the nontangential limit in (3.23) for all sufficiently large t is part of the 
conclusion of the theorem. 

Proof. We omit the proof since it is similar to the proof of an analogous theorem in [18] . □ 

Remark 3.2. When applying Theorem |3.1| to the mKdV equation in Section [4j we will identify 
( with the quotient ( := x/t and choose X = (0,1V] with N > 0. Moreover, the two critical 
points will be located at ±fco with fco := \J C/12 and we will take e := fco/2 and p := ei/48k(, so 
that t = tp 2 = ^/x 3 /(\21). The contours T and f will be as displayed in Figure [ 6 ] The idea is 
that the disks of radius e centered at ±fco provide small neighborhoods of the two critical points 
in the fc-plane, such that as t —> oo the jump matrix is close to the identity matrix everywhere 


except on the two small crosses ±fc 0 + X e (cf. assumption (3.17)). Near the critical point fc 0 , 
it is convenient to work with the variable z = (j. e . k — — £ 2 ) which centers the 

critical point at the origin and maps the cross fco + X e of radius e in the fc-plane to the cross 
X p of radius p in the z-plane. The scaling factor p is chosen so that the jump matrix takes the 
standard form (3.19) in the 2 -plane. The variable t is introduced so that t —» oo corresponds 
to the condition x 6 ft 


oo in the definition ( 1.21 of the similarity sector. 


4. Asymptotics in the similarity sector 
The goal of this section is to prove Theorem |4.1| which provides an asymptotic formula for 


the quarter plane solutions of the mKdV equation (1.1) in the similarity sector. The proof is 
essentially an application of Theorem |3.1[ 


(4.1) 


Let £ = x/t and define the variables fco = fco(C) and r = T i x U) by 


kn = \ — 


12 ’ 


t = 12 tk. 


o- 


Theorem 4.1 (Asymptotics in the similarity sector). Suppose r{k) is a function in C 11 ( 
such that 


• r{k ) = r(—k) for k € R. 

• sup feeR |r(fc)| < 1 . 

• r^\k) = 0 (k ~ 4+2j ) as \k\ —t oo, k € R, for j = 0, 1 , 2. 

Suppose h{k ) is an analytic function of k £ D 2 such that 

• h^\k) has a continuous extension to D 2 for each j = 0,1,..., 5, that is, h € C 5 (D 2 ). 

• h(k ) = h(—k) for k € D 2 . 


There exist complex constants {hi}f such that 


h^\k) = 


di 

dkl 


hi 

k 


+ pl +oyr 4 «i), 


oo, fc € dDi, j = 0 , 1 , 2 . 


Let a £ [1/2,1) and let N > 0. 


Then there exists a T > 0 such that the L 2 -RH problem (2.9) has a unique solution and the 


limit in (2.10) exists whenever r > T and 0 < x < Nt. Moreover, the function u(x,t) defined 
by ( 2 . 10 ) satisfies 


1 
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where the error term is uniform with respect to x in the given range and the function u a is 
defined by 




with 


V’(Cs) = 


-- \j 3 cos (l6ffcg ln(192ffcg) + </>(C)) 

(4.3) 

+ argT(ii/(£)) arg r(k 0 ) + i ^ j , 

Tjr s k 0 

(4.4) 

7 ^1n(l - Hko)\ 2 ), 

(4.5) 

fln(l - |r(s)| 2 ), |s|>fc 0 , 

[ ln(l - |r(fc 0 )| 2 ), |s| < k 0 , 

(4.6) 


Remark 4.2 (Quarter plane solutions). Theorem 4.1 holds for any functions r(k) and h(k) 


which satisfy the stated assumptions, regardless of whether these functions are defined via (2.3) 
or not. By assuming that r(k) and h(k) are defined in terms of uq(x) and {gj(t )}q via (2.3) and 


by combining Theorem 2.1 and Theorem 4.1 we can ensure that the function u(x,t) defined in 


(2.10) constitutes a quarter plane solution of (1.1) with initial and boundary values given by uq 
and More precisely, we have the following corollary: Suppose Uq, f/ 0 , , g - 2 satisfy (2.4) 


with n= 11 and m = 4, that t he a ssociated spectral functions satisfy the global relation (2.51, 
and that the L 2 -RH problem (2.9) has a solution for each (x,t) £ [0, oo) x [0, oo). Then the 


function u(x,t) defined by (2.10) is a well-defined classical solution of (1.1) in the quarter plane 
with initial data u(x, 0) = uq(x) and boundary values dfu(0,t) = gj(t), j = 0,1,2. Moreover, 
the asymptotics of u(x,t) in the similarity sector is given by (4.2). 


Remark 4.3 (Effect of the boundary). According to equations (4.3)-(4.6), the asymptotic 


wave shape u a (x,t) associated with a quarter plane solution u(x,t) depends on the initial and 
boundary values via the function 


r{k) = 


b(k) B(k) 


a(k) a(k)d(k )’ 


where d(k) = a(k)A(k ) — b{k)B(k). The functions {a(k), b(k)} are defined in terms of the initial 
data, whereas the functions {A(k),B(k)} are defined in terms of the boundary values. If we 
considered equation (1.1) on the line (instead of the half-line), the solution would satisfy an 
asymptotic formula analogous to (|4.2|) but with r(k) replaced with 


\k) = 


b L (k ) 
a L (k) ’ 


where {a L (k), b L (k)} are the line analogs of the half-line spectral functions {a(k), b{k)}. It 
follows that the effect of the boundary on the leading-order asymptotics in the similarity sector 
is to add the term — B{k)/{a(k)d(k )) to the value of r(k). 


Proof of Theorem \4-l\ L et N > 1 be given and let I denote the interval I = (0, N]. The jump 
matrix J defined in ( |2.8[ ) involves the exponentials where 4)(Qfc) is defined by 

$(C, k) = 8 ik 3 — 2 ikf, f £ X, k £ C. 

It follows that there are two stationary points located at the points where 'jyf = 0, i.e. at 
k = ±fc 0 . The real part of is shown in Figure [3j 

Step 1: Introduce an analytic approximation of h(k). The first step of the proof 
consists of decomposing h into an analytic part h a and a small remainder h r . 

Lemma 4.4. There exists a decomposition 

h(k) = h a (t,k) + h r (t,k), t > 0, k £ dDi, 
where the functions h a and h r have the following properties: 
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Im 



Re k 

Figure 3. Contour plot o/Re<£’(£, k) for ko = 1. Light (dark) regions correspond to positive 
(negative) values of Re$. 


(a) For each t > 0, h a (t , k ) is defined and continuous for k £ D\ and analytic for k £ D\. 

( b ) For each ( £ T and each t > 0, 

C 


\h a (t,k)\ < 


e f|Re<6(C,fc)| 


h G D \, 


i + |fc| 

where the constant C is independent of (,t,k. 

(c) The L 1 , L 2 , and L°° norms of the function h r {t , •) on dD\ are 0(t~ 3 / 2 ) os t —> oo. 

(d) The following symmetries are valid: 


h a (t,k) k)j h r (t, k) h r (t , kf 

Proof. Since h £ C 5 (D 2 ), there exist complex constants { pj }q such that 


(4.7) 


v j=o 


0 , k £ dD i, n = 0 , 1 , 2 . 


In fact, pj = , ‘ 1 '-, (0) 


for j = 0,1,... ,4. Let 


fo(k) = 


i=i 


(fc — i)i 


where {a.,}® are complex constants such that 


Mk) = l ZUPi kj + o( k ^ 

M ^E -=1 h 0 k-i+ o (k-% 


o, 

OO. 


(4.8) 


It is easy to verify that (4.8) imposes eight linear conditions on the a/s that are linearly inde¬ 
pendent; hence the coefficients aj exist and are unique. The symmetry h(k) = h(—k) implies 
that aj = (—iydj, j = 1,..., 8. Hence fo(k) = fo(—k). Letting / = h — / 0 , it follows that: 
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(i) fo(k) is a rational function of k G C with no poles in D\. 

( ii ) fo(k) coincides with h(k) to order four at 0 and to order three at oo; more precisely, 


d n f 

dk n 


(k) = 


0 (k 5 ~ n ), 0, 

0 (k~ i+2n ), k —> oo, 


k G 9Di, n = 0,1,2. 


(4.9) 


(in) f(k) = f(—k) for k G C. 

The decomposition of h(k) can now be derived as follows. The map k <—► (j> = 0(/c) defined by 

0(fc) = 8fc 3 

is a bijection dD\ -A M, so we may define a function F : K. —► C by 


F(<j>) = (k — i) 2 f(k), CeX, 


(4.10) 


The function F(</>) is C 5 for <fi y 0 and 


F (i) (0) = 


i a 

24 ¥dk 


[(k-i) 2 f(k)\, <t> G K\ {0}. 


By (4.9), F G C' 1 (K) and F"(<j)) = 0(\(j>\ 1/<3 ) as (j) —> 0, while F^(0) = 0(101 2 / 3 ) as |0| — > oo 

<oo, j = 0,1,2, 


for j = 0,1, 2. Hence 


d J F 

dty 


m 


that is, F belongs to the Sobolev space FT 2 (R). It follows that the Fourier transform F(s) 
defined by 


satisfies 


and 


H s ) = y i F^y-^d^ 

J M 

F(0) = f F(s)e i<t,s ds, 

J R 


(4.11) 


Equations (4.10) and (4.11) imply 

1 / 


(.k - i) 2 


|s 2 F(s)|| L 2 (r) < oo. 


F(s)e 8isk ds = f(k), k G dD\. 


Writing 

f(k) = f a (t, k ) + f r (t, k), t> 0, k G dD u 
where the functions f a and f r are defined by 


fa(t,k) = 
fr(t,k) = 


1 


(k-i ) 2 J _| 

1 


F(s)e 8 ' lsk ds, t > 0, Jfc G Di, 


F(s)e 8xsk da, t> 0, k G &Di, 


(fc - i ) 2 J_ c 

we infer that / a (f, •) is continuous in Di and analytic in D\. Furthermore, since |Re8ifc 3 | < 
|Re$(0, k )| for all A: G D\ and ( £ 1, we find 


< 


c 


\k~i\ 


e !IRe$(c,fe)| 


£ G X, £ > 0, A) G D\ 
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and 

— - I — 1 

\fr(t,k)\ < | fe ^ .|2 J s 2 \F(s)\s~ 2 ds < 1 ^—T ^\\s 2 F(s)\\ L 2 {m ^j J s~ i ds 

- i + \ k \2 rV2 ’ *> °> k£dD 1 . 

Hence the L 1 , L 2 , and L°° norms of f r on dD\ are 0(t~ 3 / 2 ). The symmetry f(k) = f(—k ) 
implies F(<)>) = F(—(f>). Thus the Fourier transform F(s ) is real valued, which leads to the 
symmetries in ( |4.7| ). Letting 

h a (t,k) = f 0 (k) + f a {t,k), t> 0, k £ D ± , 
h r (t , k) = f r {t , k), t > 0, k £ 9Z?i, 

we find a decomposition of h with the properties listed in the statement of the lemma. 

V 


Step 2: Deform. Lemma |4.4| implies 

fc«(f>-)e** (c, ° £E 2 {D 1 )nE°°(D 1 ) 

for each £ £ X and each t > 0, where E°°(D i) den otes the space of bounded analytic functions 
in D\. Thus, M satisfies the L 2 -RH problem (2.91 iff the function M h (x,t, k) defined by 


M h {x 1 t 1 k) = < 


satisfies the L 2 -RH problem 


M(x,t , k) 

M(x,t , k) 
k), 


1 ° 

-h a {t,k)e t9 ^ 1 J 

.0 1 


k £ £>i, 

k £ D 4 ., 
otherwise, 


(4.12) 


,•) e / + f; 2 (c\e), 

M+(x,t,k) = M k (x,t,k)J h (x,t,k) for a.e. k £ E, 


(4.13) 


where 


t, fc) = < 


1 ° 

K h r (t,k)e t ^’ k '> l)' 

1 — r(k)e~ t ^t k ) s 

. r(fc)e t$ ^ ,fe - ) 1 —|r(fc)| 2 ^ 

1 h r (t,k)e- t *tt’ k A 

,0 1 J ’ 


/c £ dDi, 

k £ R, 
fc G 9I?4. 


(4.14) 


Part (c) of Lemma 4.4 implies that the L 1 , L 2 , and L°° norms of .J h (x. t,-) — I are 0(f 3 / 2 ) on 
dD 1 UdD 4 . 


Step 3: Conjugate. In order to apply the steepest descent result of Theorem |3.1[ we 
need to transform the RH problem in such a way that the jump matrix has decay everywhere 
as t —» 00 except near the two stationary points. This can be achieved by performing an 
appropriate triangular factorization of the jump matrix followed by a contour deformation [lOl . 
For k £ (—fco,fco), ft eas y to achieve an appropriate factorization. We next show that an 
appropriate factorization can be achieved also for k £ (—00, —ko) U (A'q, oc) by conjugating the 
RH problem ( 4.13| ). 

Let 

'S(Ck)- 1 0 \ 

0 S(C,k)J’ 


A(C, k) = 
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where 

S(( ,/e) = +/^ )ln C 1 H r(s ) |2 )^ ; k £ C\R. 

The function S satisfies the following jump condition across the real axis: 


* + (c,fc) = |fc|>A:o ’ 

\S-(C,k), \k\<ko, 


k e 


Moreover, the symmetry r(/e) = r(—k) implies that 

5((,k) = 6(<;,-k) = 5((,-k)- 1 . 
It follows that A obeys the symmetries 

A(C) k) = A(C, —k) = A(C, —k)^ 1 . 

Lemma 4.5. The 2 x 2-matrix valued function A(£, k) satisfies 


A(C,-),A (C,-)" 1 e/ + s s 


n£° 


/or eac/i £ £ I. 


(4.15) 


(4.16) 


(4.17) 

(4.18) 


Proof. First note that 


where v{(f) is given by (4.5), 




X((,k) = f 

Jr 


e x(C,fc) ; 

(4.19) 

ds 

(4.20) 

s) s-k ' 


and the function s) is defined by (4.6). Since •) £ H 1 ( R), we have x(Ci ') £ E° 
for each / £ Z; in fact, see Lemma 23.3 in [2], 

sup sup |x(C,z)l < sup||i/)(C,-)II hi(R) < °°- 

CeizeC\R cex 


(4.21) 


Hence S((, •),£(£, -) -1 £ E °°(C \ R) for each £ £ I. 

It remains to show that <5(£, ')> <5 (Cj ’) _1 £ 1 + ,E 2 (C\R). Since ip(Cr) £ L 2 (R), x(Cc) £ 
E 2 (C \ R) for each / £ I. Let C n be the image of the circle |w| = 1 — i under the linear 
fractional transformation z K > —i z f l , which maps the unit disk onto the upper half plane. 
Then sup„ ||x(C> ')l|i 2 (C„) < Using the inequality 

\e w — 1| < |u>| max(l, e Rew ), w £ C, (4.22) 


we obtain 

sup ||A(C, •) - l||z,2 (Cn) < sup 


k — k 0 
k + k 0 


L 2 (C n ) 


+ sup \\e x - l||z, 2 (Cii) 


< C sup ||e x ||ioo (Cri) + sup ||xe |Rexl ||i 2 ( c> ri) < oo. 

n n 

It follows that S((, •) £ 1 + E 2 (C \ R). A similar argument gives <5(£, -) _1 £ 1 + E 2 (C \ R). V 


By Lemma [4.5[ M h (x,t,k ) satisfies the L 2 -RH problem 

M{x,t,k) = M h {x,t,k)A((,k) 

satisfies the L 2 -RH problem 

M(x,t ,■) £ / + E 2 (C \ X), 

M + (x,t,k) = M-(x,t,k)J(x,t,k) for a.e. k £ X, 


4.13) iff the function M defined by 


(4.23) 
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k G 8D i, 


where 

J[x, t, k) = Al 1 (^, k)J h (x , t, k) A + (C, k) 

' 1 o\ 

-MC, fc)MC, 

,r(fc)e t$(C ’ fc) f^|(l - k(fc)| 2 ) , 

k G 9Z?4. 




■MC.fc) 

<5 + (C,fe) 

^-(c,fc)V(c.fc) rrfeV ** (c,fc) 


k G R, (4.24) 


A <S(£ ,fe) 2 /i r (^fc)e- t$ ^’ fc ) 

,0 1 


In view of the jump (4.16) of <5(£, fc), this gives, for k G R, 


1 - |r(t)| 2 
Ix.CG r ( fc ) _ 

+ l-|r(fc)|= e 


1 


, \k\ > ko, 

\k\ < k 0 . 


1 —5((,k) 2 r(k)e 

\8{C,,k)~ 2 r{k)e t ^^' k ' ) 1 — |r(fc)| 2 J 

The upshot of the above conjugation is that we can now factorize the jump matrix as follows: 


J = 


B u 1 B t , |fc| > fc 0 , k G R, 
b~[ 1 b u |fc| < k 0 , k G R, 


(4.25) 


where 

. 1 

y 6 + ((,k) 2 n(k)e t ^U, k ) lj ’ 

1 0\ 

K —S((, fc) _ 2 r 3 (/c)e t$(< ’ ,fc) 1 J’ 

and the functions {fj(k)}f are defined by 

r(k) 


bu — 


1 — <$(£, fc) 2 r 2 (fc)e 


0 


1 


bi = 


_/l <5_ (C, fc) 2 r4(fc)e _t$( '»’ /c) 

“ _ lo i 


(4.26) 


ri(fc) = 


1 — r(k)r(k) 
r 3 (k) = r(fc), r 4 (fe) = 


r 2 (fc) = r(k), 


r(k) 


1 — r(k)r(k ) 


Our next goal is to deform the contour near R. However, we first need to introduce analytic 
approximations of {r ? (fc)}j. 

Step 4: Introduce analytic approximations of {rj(k)}f. We introduce open sets Uj = 
Uj{ £), j = 1,... ,4, as in Figure [4] so that 

{k | Re <F(C, k)<0} = U 1 U U 3 , {k | Re $(£, k) > 0} = U 2 U U A . (4.27) 


The following lemma describes how to decompose rj, j = 1,... ,4, into an analytic part ry a 
and a small remainder rj iT . A proof can be found in [lftj . 

Lemma 4.6. There exist decompositions 

rj, a {x,t,k)+rj tr (x,t,k), j = 1,4, |fc| > k 0 , k G R, 


rj{k) = 


rj, a {x,t,k) +r jt r(x,t,k), j = 2,3, |fe| < k 0 , fcG 


where the functions {‘>'j,aji"j,r}j—i have the following properties: 

(a) For each £ G T and each t > 0, rj^ a (x, t, k ) «s defined and continuous for k G Uj and analytic 
for k G Uj, j = 1,... ,4. 
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Ua 


Figure 4. The domains {Uj}f in the complex k-plane. The dashed curves are the curves on 
which Re <£> = 0. 



Figure 5. The contour T and the open sets {Vj} 


(&) The functions {rj, a }i satisfy, for each I\ > 0, 

\r jt a(x,t,k) - rj{ka)\ < C K \k - fc 0 |e5 |Re<I,(c ’ fc)l , 

k&Uj, \k\<K, (el, t > 0, j = 1,... ,4, 

where the constant Ck is independent of C,t,k but may depend on K. 

(c) The functions and r^ a satisfy 


\rj, a {x,t,k)\ < 


C 


1 + |fc| 


e f|Re$(C,fc)| 


k € Uj, ( € I, t > 0, j = 1,4, 


(4.28) 


(4.29) 


where the constant C is independent of(,t,k. 

( d ) The L 1 ,L 2 , and L°° norms on (— oo, —ko)U(ko, oo) of the functions ri tT {x,t, ■) andr^ r {x,t, ■) 
are 0(t~ 3 / 2 ) as t -4 oo uniformly with respect to £ € I. 
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(e) The L 1 ,E 2 , and L°° norms on (—ko,ko) of the functions r 2 , r {x,t,-) and r 3ir .(a;, f, ■) are 
0(t~ 3 / 2 ) as t —► oo uniformly with respect to £ G I. 

(/) The following symmetries are valid: 


= r iAC,,t,-k), r j>r (C,t,k) = r jiT (C,t, -k), 


j = !,•••, 4. 


(4.30) 


Step 5: Deform again. Let T be the contour consisting of E together with the four 
half-lines 


{ko + ue ^ 1 4 | — V 2 k 0 < m < oo }, { — &o + ■ 


± dl7T 
A 


— V 2 ko < u < oo}, 


Bu B u , r B u , a , 


oriented as in Figure[5j Let {V}}® be the open sets shown in Figure|5] Write 

Bi = Bi r Bi a , b u b ur b ua , bi = bi r bi a ? 

where {B^ a , 6 U|Q , b Ua , B Uta } and {Bi >r , b u>r , bi >r , B u>r } denote the matrices {B t , b u , b u B u } with 
{rj(k)}f replaced with {rj ta {k)}f and {rj, r {k)}f, respectively. The estimates (4.281 and (4.29) 
imply that 

r B l:a (x,t, •), B-'(x,t, •) g / + (e 2 n £°°)(Hi), 
b Uta (x,t,-),b~^ a (x,t,-) g / + {e 2 n E°°)(y 3 ), 
bi,a(x, *. •)> *,■)€/ + (£ 2 n £°°)(h 4 ), 

^B u ,a(x,t,-),B-i(x,t,-) g / + (^nF®)^), 


for each (£l and each t > 0. Hence M satisfies the L 2 -RH problem (4.231 iff m{x, t, k ) defined 
by 


i(x, t, fc) = < 


satisfies the L 2 -RH problem 

f mp, t, •) G / + F? 2 (C \ T), 

[ m+p, t, fc) = m~{x, t, k)v(x, t, k) for a.e. /c G F, 


M(x,t, k)B i;a (x,t, k)~ l , 

k G P, 

M(x,t, k)b Ui a{x,t, k)~ l , 

* G p, 

M(x,t, k)biAx,t, ky 1 , 

k G V4, 

M(x, t, k)B u y x, t, ky 1 , 

fc G H 6 , 

M(x, t, k), 

elsewhere, 


(4.31) 


where, in view of (4.24) and (4.25), the jump matrix v is given by 

1 


v = 


Bl n — 


O' 


l,a \S(C,k) 2 r i,a{x, t, k)e t<s? ^ ,k ^ 1 j 

( 1 -<KC, k) 2 r 2 , a {x, t, fc)e _t$(c ’ fc)N \ 

j 

1 ° I 

v -<5(C,fc)“ 2 r 3 , a (a;,t,fc)e t<I,(c ’ fe) 1J 
h- <5(C> k) 2 r 4 , a (x,t, fc)e _t$(c ’ fe)N \ 


bu,a 


bl.a = 


3 u,r a l,r j 


b, r b u , r , 

1 O' 

<5(C, k)~ 2 h r (t, /c)e t$ ^ ,fc ^ iy 
1 (J(£, k) 2 h r (t, 

0 1 


k G V\ n V2, 
fe g H 2 n H 3 , 
/j g V 4 n V5, 

kp h 6 . 

k G p n h 6 , 
fc g p n V4, 

k G 91>i, 
fc G dD 4 . 


(4.32) 
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Figure 6. 


The contours T and T satisfy the properties (n)-(r4) of Theorem 


3.1 


From the symmetries (4.71, (4.17), and (4.30), we infer that v satisfies 

v(x, t, k) = v(x, t, — k), k€ r. 


(4.33) 


Step 6: Apply Theorem |3.1[ We verify that Theorem |3.1| can be applied to the interval 
X = (0, N], the contour T, and the jump matrix v with 


k 0 

e ~ T’ 


/ 52$ ,— 

= = e v 48fc o, 


9(C) = e 


if 

- (J 7T l Jl 


1 £ Tp(C,s)ds 


~r{k 0 )e 


2iv{C,) ln(2 v / 48fco /2 ) 


r = tp 2 = 12 /cq^, 

KO = ln (! - Hko)\ 2 ), 


KO z) = $ ( C, k 0 - -z ) = —16 */cq + 


IZ 0 

12p' 


The contours F and f are shown in Figure [6] The conditions (n)-(r4) are clearly satisfied. 
Since the contour is independent of (, a scaling argument shows that ll^f IIb(l 2 (T)) ' s 


independent of In particular, Sp is uniformly bounded on X 2 (r). Equation (|3. 14 ) follows 
from (4.32) and the estimates in Lemmas 4.4 and 4.6 Clearly deti; = 1. The symmetry 
condition (3.16) follows from (4.33). 

We next verify (3.17). Let w = v — I and let R denote the contour obtained from T by 
removing the crosses ±ko + X e . By parts (d) and (e) of Lemma 4.6 the L 1 and L°° norms of 


are Oft 2 ) < 0(er x ) on R uniformly with respect to Q G X. By part (c) of Lemma 


4.4 


the 

L 1 and L°° norms of w are 0(f _2 ) < 0(er _1 ) on dDi U dD 4 uniformly with respect to £ G X. 
Let 7 denote the part of V that belongs to the line ko + Re 1 ?, i.e. 


7 = < fco + we 4 


u G 


-V2k 0 ,- k 4 


k 0 

2 ’°° 


Let k = ko + ue 4 . Then 


Re $(0 k) = — 4w 2 (6fco + V%u) < —16u 2 ko for — V2k 0 < u < 00 . 


(4.34) 


Together with (4.29) this yields 


„ pi 2r. 12tu 2 ko 

\r ha (x,t,k)e t,5> ^\ < C7 e -¥lH»*(C,*)l < Ce~ 12t ^ k ° < CtU ^ - 

tu z ko 


< 


C 


< 


Ck 2 


k 0 

— < u < 00. 
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Similarly, by ( |4.28 l, 


|r 3ia (x,f,fc)e t$(c . fc )| < Ce ~^ |Re<I,(C ’ fc)l < 


Ck 2 k 0 

—S—, -V 2 kn<u<- — . 


Hence 


Ck 2 


klkbT) ^ 


\w\\ L oo M =°{t ). 


-V-!■'(■) 


u 2 du = O(k 0 r ), 


This shows that the estimates in (3.17) hold also on 7. Since similar arguments apply to the 
remaining parts of T', this verifies (3.17). 


Equation (4.32) implies that (3.18) and (3.19) are satisfied with 

= <5(C ,k)~ 2 r 3ta {x,t,k)z 2w ^\ 
R 2 {C,t,z) = <5(C, k) 2 T4 t a{x, t, k)z~ 2lV ^\ 
#3(C>M) = <00 k)- 2 r lt a{x,t,k)z 2w ^\ 
^R A (C,t,z) = S(C,k) 2 r 2 ,a{x,t,k)z~ 2w( ' <: '>, 

where k and z are related by k = ko — —. 


It is clear that (j> satisfies (3.20) and (3.21c). The following estimate proves (3.21b): 

U|2 


Re 0(00 = 


1 ± 


6^2 p 


> 


4 ’ 


&x p ux p , (el, 


where the plus and minus signs are valid for z G X 2 and z G X p respectively. The proof of 
(3.21a) is similar. Since |</(0| = |r(fco)|, we have sup^ e:r |(?(0| < land^(C) = — ^ ln(l — |g(C)| 2 )- 
Finally, we show that given any a G [1/2,1), there exists an L > 0 such that the inequalities 
(3.22) hold. Let k = ko — —. Using the expression (4.19) for 6(C,k), we may write 

=e- 2x(c ’ fe) r 3 , 0 (x,t,A:)e 2il/(c)ln( (- fcc> - fe )' /3SSo ), z € X,. 

Thus, 

Ri(C,t, 0) = e-^- ( «- feo V 3ja (x,t,fco)e 2w(c)(ln(2 ^ fe ° 3/2)+l7r) . 

Now r 3 i a(x,t,k 0 ) = r(fco) by ( |4.28[ ), and 

X-((,ko ) = -x— f - 7 -MO 

2 7r *jR s — ko 

by Sokhotski-Plemelj. Hence I?i(C,f, 0) = (/(C)- Similarly, we find 


-MOM) — 


9(0 


1-|9(0I 2 ’ 


MOM) = 


9(0 


1-|9(0I 2 ’ 


MOM) =9(0- 


We establish the estimate (3.22) in the case of z G X p ; the other cases are similar. Note that 
z G X p corresponds to k G k 0 + XJ. 

We have, for z G X p , 


IM0M) - q{ 01 < 

o -2 x _(C,fe 0 ) 


e -2 x((,k) _ e -2x-(C,fco) 


|r 3 , a (x,t, fc) - r 3 (fc 0 )l 


r 3 } a{x, t , fc)e 2MC) in((-fco-fc)V®S5) 
2ii/'(£) ln((—fco - fc)\/48fco) 


e -2 x _(C,feo ) r3 ( fco ) 


a 2 ii 7 (^)(ln((fc 0 +fe)\/ 48 fco)+iir) _ e 2 ii 7 (C)(ln( 2 v / 48 fcg^ 2 )+i 7 r) 


The functions e 2 *M) ln (( fc ° fe ) v/48fe °) and e 2 x-(C,fco) are uniformly bounded with respect to 
k G ko + X 3 and C € 1. Moreover, employing the estimate 


|Re</(C,ue")| = 


-- 1 - 


072 pj 


2 v 2 

< - 3 -, -P<v<p, 
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we see that equation (4.28) yields 


|r3,„(*,i, A;) - r 3 (fc 0 )| < C\k - fc 0 |e* |Re$(c ’ fc)l = 

P 


P 


G Xf. 


Thus, 


e -2x(C,fe) _ e -2x_(C ,ko) 


ez LiT 

+ C—e 6 
P 


\Ri((,t,z)-q(0\<Ce-r- 

+ C \i- e - 2i ^ ln ^ ) \, (el, t > o, zex{. 

The estimate in (3.22) for z G X{' now follows from the following lemma. 

Lemma 4.7. TVie following inequalities are valid for all C G X and all k G ko + XJ: 

< (7|fc — fe 0 |(l + | In |fc — fc 0 ||), 


e -2 X (C,fc) _ e -2x-(C,feo) 


1 — e 


-2iK01n(^l 


<C'feo 1 |fc-*o|, 


where the constant C is independent of C and k. 

Proof. We first prove that 

|X(C)A:) - X-(C,A: 0 )| < C\k-k 0 \(l + |ln|fc - fc 0 ||), C G X, k G fc 0 + XJ. 
Integration by parts in the definition ( |4.20| ) of x gives 

X(C) k) = ( J + ^ ln(s - fc)dln(l - |r(s)| 2 ). 


Hence 


lx(Ci k) — X-(C) & 0 )| = 


27T 


—fco /»oo 
k 0 

p — ko poo 

+ / 

— OO -7 fco 


In 


This gives 


I x(C,k) - x_(C,fc 0 )| ^ C (yJ 


r-k 0 rco \ 

, s — k 

/ + / 

in 

J—OO j /cq / 

s - k 0 


^-A)<;i»(i - k(0l 2 ) 

ds 


< C 


/ 2/cq 


In 1 + 


k — kr 


du 


C 


l + |s| 2 

lull-Git? 


The change of variables v = \k — k 0 \/u yields 
k - k 0 


In 1 - 


1 + (u- k 0 ) 2 J 0 

0 \/u yields 

) 7 POO 

= \ k - fc °l J o |ln(l+t ; eT)| 


du 


1 + u 2 

dv 


v 2 + \k — ko\ 2 


Since 


|ln(l + ve'*) \ < C 


v , 0 < v < 2, 

lnii, 2 < v < oo, 


we conclude that 


In 1 - 


k - kr 




vdv 


lnw 


v 2 + | k - fc 0 | 


dv 


< C\k - /sold In I k - foil + C), k€k 0 + X£. 


A similar argument applies to the first term on the right-hand side of (4.40). This proves 
Using the inequality (4.22) together with (4.21) and (4.38), we estimate 
|g-2 X (C,fc) _ g—2x— (C.*=o) I < |g-2x-(C,fco)||g-2[ X (C,fc)-X-(C,fco)] _ ]_| 

< C\x(C, k) - X- (C, A: 0 ) |e 2|Re ^’ k )-x-(CM))\ 

< C\k — A;o|(l + | In |As — fc 0 ||), C G 2, A: G 1 + X|. 


(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 


(4.38). 
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This proves (4.36 1 . 


By (4.22), 


1 — e 


- 2 MC)ln(%±A) 


< 


2z/(C) In 


< C 


In 1 + 


k 0 + k 
2 k 0 
k — k 0 

2ki 


o 


which proves (4.37). 


JRe (2^(0 ln(*a+*))| 


< Ck 0 1 |fc — fc 0 |, keko + X'. 


3’ 


V 


Step 7: Find asymptotics. As a consequence of Theorem 3.1 the I/-R H problem (4.31) 


for m has a unique solution and the limit in (3.23) exists for all sufficiently large r. Since the 
RH problems for M and m are equivalent, the /A-Ri 1 problem (2.9) for M also has a unique 
solution for all sufficiently large r. Moreover, we conclude that the limit (2.10) defining u(x,t) 
exists whenever r is large enough and 

z z 

u(x,t) = — 2 i lim (kM(x,t,k)) 12 = — 2 * lim (km(x,t,k)) 12 - 


k—> 00 


k—> 00 


Equation (3.23) of Theorem 3.1 then yields 


u(x,t ) = — 


4eRe/?(£,£) 


\ Re/3(C,£) / _ i+Q 

OUt 2 )=- ;’ + qUt 2 

V ’ VZtko y 


where (3 (C, f) is defined by (3.24). We have 

Re/3(C ,t) = \J 12 (C) cos - arg 5 (C)+ argT(w(C)) + 16£fco- 12 (C) In 


where 


Thus, 


arg q{() = 2p(C) ln(2/c 0 \/48fco) + argr(fc 0 ) - 


^((,s)ds 
S - fc 0 


u(x,£) = 


'K0 

3tA : 0 


cos(16£fcQ — 12 (C) ln(192£fcg) + 0(C)) + 0(er 2 ), 


where 0(C) is defined by (4.4). This proves (4.2 1 and completes the proof of Theorem 4.1 


□ 


5. A MODEL RH PROBLEM 


The remainder of this paper is devoted to deriving the long-time asymptotics for the mKdV 
equation in the self-similar sector ( |1.4| ) . The goal of the present section is to prove Theorem 
PI which expresses the large z behavior of the solution of a model RH problem (see equation 
(5.2)) in terms of the solution of the Painleve II equation. The result of Theorem 5.1 will be 
important for the proof of Theorem |6.1| 

Let Z denote the contour Z = Z\ U Z 2 U Z 3 , where the line segments 


Zi = {1 + re " | 0 <r<oo}u{ —1 + re 5 " | 0 < r < 00 }, 

Z 2 = {—1 + re - ^ |0<r<oo}u{l + re~ t | 0 < r < 00 }, 

Z 3 = {r | - 1 < r < 1} (5.1) 

are oriented as in Figure [7j Given a positive number zq > 0, let ZqZ = {zqZ € C | z £ Z} denote 
the contour Z scaled by the factor zq. It turns out that the long time asymptotics in the self- 
similar sector is related to the solution m z of the following family of RH problems parametrized 
by s G *K, y G M, Zq > 0: 


m z (s,y,-,z 0 ) G I + E 2 (C\z 0 Z), 

m z (s, y, z, zq) = m z (s,y,z,z 0 )v z (s,y,z,z 0 ) for a.e. z G z 0 Z , 


(5.2) 
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Figure 7. The contour Z. 


where the jump matrix v z (s, y, z, ~o) is defined by 

1 °^l 

se 2i(yz +±#) 

z , . ./I - se - 2i (^+¥) N 

V ( S,y,z,z 0 ) = , ( 

1 se~ 2i ( yz+i t“) N 

0 1 


1 


O' 


e 2 i(yz+^) i 


z g 0 O ^i, 
z G ^0-^2, 
2 € ZqZ 3 . 


(5.3) 


Note that —> I exponentially fast as 2 -> oo, z G ZqZ. Moreover, v z obeys the symmetry 


v (s, y, z, z 0 ) = a 3 v z {s, y , - 2 , z 0 )a 3 . 

The following theorem expresses the asymptotics of m z in terms of the solution u p of the 
Painleve II equation 

(5.4) 


P \3 


Uyy = yu + 2 (u ) 

The orem 5.1. Let s G iR be such that |s| < 1. Let y G R and zq > 0. Then the L 2 -RHproblem 
{'5.2) has a unique solution m z (s,y, z, zq) which satisfies 


m z (s , y, z, zq) = I + ( P 


\y;s,0,-s) 


2 z \u p (y,s, 0 ,-s) 


Ol-U 


oo, 


(5.5) 


where u p {-\ s, 0, — s) denotes the solution of the Painleve II equation (5.f) corresponding to 
(si, S 2 , S 3 ) = (s, 0, —s) according to (A. 6), and the error term is uniform with respect to zq > 0, 
arg z G [0, 2tt], and y in compact subsets of R. Moreover, if K is a compact subset o/R, then 


sup sup sup \m z (s,y, z, zq)\ < 00. 

z 0 >OyGK z&C\z 0 Z 


(5.6) 


Proof. Uniqueness follows since de t v = 1. Moreover, the RH problem (5.2 1 for rn z can be 
transformed into the RH problem ( |A.3[ ) for m p with jump data specified by (si,S 2 ,S 3 ) = 
(s, 0, — s). Indeed, let {Oj}§ be as in Figure [ 8 j For each y G R and z 3 > 0, we have 


e-'-' rJ G (^nU 00 )^), e 


a 2i(yz+(i?2 

Hence the function m p (y, z) defined by 

'm z (s,y,z,z 0 ) 


-2 ilyz+^r) G (E2 n E°°)(n 2 ). 


z G Ho, 


m p {y, z) = < 


m z (s, y, z, zq) 


1 O' 

'l se -2i(^+ i #) N 


m z (s,y,z,z 0 )[ se2 .- yz+ ^ ) i ), zeSh, 


0 


1 


£ G ft 2 , 
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Figure 8. The sets Tlj, j = 0,1, 2. 


satisfies the RH problem ( A.3|) with jump data specified by (si,S 2 >S 3 ) = (s, 0, — s) if and only 


if satisfies the RH problem 


Since 


S £ 


i nter sect R. for S = ( s , 0, —s), see |15j and p. 51 of 
(5.6 1 follows from the properties of m p . 


and |s| < 1, the set Y$ in (A.5) does not 
[161. Thus ( |5.5[ ) follows from (|A.5|). Equation 

□ 


6. Another nonlinear steepest descent theorem 


In the self-similar sector S defined in (1.4), the asymptotics of the mKdV equation is charac¬ 
terized by the solution of a RH problem whose jump matrix has two critical points that merge at 
the origin as t —> oo. In this section, we prove Theorem |6 .1 1 which provides an implementation 
of the nonlinear steepest descent method suitable for analyzing this type of situation. Before 
stating the theorem, we describe the set-up and make a number of assumptions. 

As in Section |4j we assume the variables £ and ko are given by 

C = x/t, k 0 = \J C/12. 


Note that 0 < tk g < C in S. Hence, as (x, t) 
approach 0 at least as fast as t -1 / 3 : 


oo within S, the two stationary points ±fcg 


0 < k 0 < Ct 3 / 3 for (x, t) £ S. 


( 6 . 1 ) 


For each point (x, t) £ S, we let Y — Y ( x , t) = k$Z denote the contour Z defined in (5.1) scaled 
by the factor fc 0 ( see Figure |9|), i.e. Y = Yi U Y 2 U Y 3 , where 


Y\ = {kg + ue T | 0 < u < 00 } U {—ko + ue 6 | 0 < u < 00 }, 


^2 = {-k 0 


0 < u < 00 } U {ko + ue 's' | 0 < u < 00 }, 


Y 3 = {u - k 0 < u < k 0 }. 


( 6 . 2 ) 


Let r = r(x, t) be a family of oriented contours parametrized by (x, t) £ S and let f = FU {|fc| = 
1} denote the union of T with the unit circle oriented counterclockwise. Assume that, for each 
point (x, t) £ S: 

(ri) F and r are Carleson jump contours up to reorientation of a subcontour. 

(r2) F contains Y as an oriented subcontour. 

(r3) F is invariant as a set under the map k 1 —> —k. Moreover, the orientation of T is such that 
if k traverses T in the positive direction, then —k traverses T in the negative direction. 
(r4) The point 00 in the Riemann sphere can be approached nontangentially with respect to 

r. 
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Figure 9. The contour Y = Y\ U 1 r 2 U Y 3 . 


Assume the operator <Sp defined in (3.12) satisfies sup^^gg ll^rll/ 3 (i 2 (r)) < 00 • 

Consider the following family of /A-R11 problems parametrized by the two parameters ( x, t) G 


\m(x,t,-) G / + F 2 (C\r), 

1 m+(x, t, k) = m-(x, t, k)v(x, t, k ) for a.e. k G F, 
where the jump matrix v(x , t, k) satisfies 

w{x,t,-) := v(x, t, ■) — i € L x (r) m 2 (r) nL°°(r), (x,t) gS, 

detu(a;, t, •) = 1 a.e. on f, (x,t)£.S, 


and 


v(x, t, k) = a 3 v(x, f, — k)a 3 , (i,() G S, k G F. 


(6.3) 

(6.4) 

(6.5) 

( 6 . 6 ) 


Suppose 


t —» 00, (x, t) G S, 


(6.7) 


I IKM,-)||iP(r\y) = 0(t 3), 

\lKM,-)IU~(rvn = 0(t"3), 

uniformly with respect to x. Moreover, suppose that for k G Y the jump matrix v has the form 

0 \ 

Ri(x,t, fc)e t$(l »’ fc) 1 J' 

1 R 2 (x,t,k)e~ t ^^’ k ^\ 

0 1 )' 

1 -iZ 4 (s,t,A;)e _t * (c ’ fc) \ ( 1 


v{x, t , k) = < 


k G Yi, 

k G y 2 , (6.8) 


O' 


,0 


1 


R 3 0M,fc)e t$(C,fe) 1. 


fc G I 3 , 


where: 


The phase <h(£, k) is given by 

<&(£, k) = — 24:ik 3 k + 8 ik 3 . 


(6.9) 


• There exist constants s G iR, |s| < 1, and L > 0 such that the functions {Rj(x,t,k)}f 
satisfy the inequalities 


\R 1 (x,t,k)-s\ < L\k\e^ Re ^’ k ^, k gYi, 
\R 2 (x,t,k) + s\ < L|fc|e3l Re<I> «’ fe )l, k G Y 2 , 
\R 3 (x,t, k) - s\ < L\k\, k£Y 3 , 

}Ri{x,t,k) + s\ < L\k\, k£Y 3 , 

for all (x, t) G S with t sufficiently large. 


( 6 . 10 ) 
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Theorem 6.1 (Nonlinear steepest descent). Under the above assumptions, the L 2 -R.H problem 
(6.3) has a unique solution for all sufficiently large t and this solution satisfies 

,p ( -X 

M 


lim (km(x,t,k))i 2 = 


U ' (3i)!/3 ’ s 


,0,-s) 


k—foo 


2(3 t)i 


+ 0 (t 3 ), t —> 00 , (x,t) £ §, 


( 6 . 11 ) 


where the error term is uniform with respect to x and u p (-; s , 0, —s) denotes the solution of the 
Painleve II equation (5-4) corresponding to (s\,S 2 ,s^) = (s,0, — s) according to (A. 6 ). 

Remark 6.2. The conclusion of Theorem |6.1| can be stated more explicitly as follows: There ex¬ 
ist constants T > 0 and K > 0 such that the solution m(x, t, k ) of (6.3) exists, the nontangential 
limit lim^_ ) . 00 (fcTO(C, t, fc))i 2 exists, and the inequality 


lim (km(C,t,k))i 2 -— 

k—f 00 


,P (r3ip/ 3;S ’®’ s ) 


2(3t)s 


< Kt~* 


holds for all (x, t) £ S such that t > T. 


Proof. Sin ce det u = 1, the solution of (6.3) is unique if it exists. Let s £ *K, |s| < 1, be the 
constant in (6.101. We henceforth suppose t is so large that k 0 < 1/2. Let m z (s, y, z, z 0 ) denote 
the solution of the RH problem (5.2) whose existence is ascertained by Theorem 5.1 We note 
that if 

— x . 1 , 

V = / 0 ,M/ 3 > z={Zt)ak, 


(3 ty/z’ 


then 


Hence we define m^x, t,k) by 


2 i[yz 


4 ^ 


=mc,k). 


—X 


m 0 (x,t,k) =m s, —— x ,(3t)3fc, (3t)3/c 0 , \k\ < 1. 

(3f)3 


( 6 . 12 ) 


(6.13) 


Since mo(x,t,-) G (E 2 AE°°)({\k\ < 1}\T), it follows that m satisfies the L 2 -RH problem (6.3) 
iff the function rh(x, t, k) defined by 

, , , t \ j m(x, t, k)nio{x, t, fc) _1 , |fe| < 1, 

m(x, t,k) = < 

\m(x,t,k), otherwise, 

satisfies the L 2 -RH problem 

(rh(x,t,-) £ I + E 2 (C \ f), 

[ rh+(x, t, k) = rh-(x, t, k)v(x, t, k ) for a.e. k £ F, 
where the jump matrix v is given by 

! mo-(x, t, k)v(x, t, k)nriQ + (x, t, k)^ 1 , fcern{|fc|<l}, 
mo(x,t,k)~ 1 , |fe| = 1, 

v(x,t,k), otherwise. 

Let us write Y\ = Y p U Yf ~, where Y p and Yf~ denote the parts of Y\ in the right and left 
half-planes, respectively. Similarly, we write Y-j = Y^~ U Yf ~, see Figure 


10 


Claim 1. The function w = v — I satisfies 

( 0(|fc|e _4 * |fc_/c ° |3 ), k£Y+UY+, 

w(x, t, k) = < 0(|fc|e _4t l fc+fc °l 3 ), k£YfUYf~, t —» 00 , (x,t) £ S, 

lo(|fc|), k£ Y 3 , 

where the error term is uniform with respect to x and k in the given ranges. 


(6.14) 
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Figure 10. Y+ and Yj denote the parts ofYj, j = 1,2, in the right and left half-planes, 
respectively. 


Proof of Claim [7J We give the proof for k £ Y+ U F 2 + U Y 3 \ the proof for k £ Yf~ U Y 2 ~ is 
similar. 

The definition (16.91) of $ implies 


Hence 


Re$(C ,k) = T(8r 3 + 12\/3k 0 r 2 ), 

j Re $(£, k) < — 8|fc — ko\ 3 , 

1 Re$((), k) > 8|fc — fc 0 | 3 , 


k = k 0 + re* 1 '® , r > 0. 
k £ Y+, 


k £ > 


(x, t) £ S. 


(6.15) 


If A; € K| + and |fc| > 1, then only the (21) entry of w is nonzero and equations (6.101 and (6.15) 
yield 

\w 2 i{x,t,k)\ = |i?i (x, t, fc)|e _t|Re$( ^ ,fe)l < (T|fc| e 5 |Re$(c ’ fe)l + |s|) e - t|Ito * (c ’ fc)l 
< (L\k\ + |s|)e- 4 ‘ |fc - fe ° |3 < C|fc|e- 4 ‘ |fe - fe ° |3 , (x,t)£ S, k£Y 1 , \k\ > 1. (6.16) 


A similar argument applies when k £ 1 2 + . This proves (6.14) for k £ Y+ U 1 2 + with \k\ > 1. 
It remains to consider the case \k\ < 1. In this case, we have 

w(x, t, k) = nio-(x, t, k)u(x, t, k)mo+{x, t, k)^ 1 , k£Y , |fc| < 1, 

where 


u(x,t,k) := v(x, t , k) — v z ( s,- r , (3 1) 3 k, (3t) 3 ko ). 

(3 t)s 


In view of equation (5.6), mo+(x,t,k) and mo_(x,t,k) remain uniformly bounded for (x,t) £ S 
and k £ Y. Therefore, it is enough to prove that 


u(x, t, k) = 


_ f0(|fc|e _4< l fc_fe °l 3 ), k£Y+\JY+, 
|o(|fc|), k£Y 3 , 


t —> 00, (x,t) £ S, |fc| < 1, (6-17) 


uniformly with respect to x and k. Equation (6.17) follows from the assumptions (6.8M6.10) 
on v. Indeed, for k £ , only the (21) entry of u(x,t,k) is nonzero and equations (6.10) and 

(6.15) yield 

3 


\(u(x,t,k)) 2 i\ = \Ri(x,t,k) - < L|fc|e- 4t l fc - fc °l , 

(x,t) £ S, k £ Y+, \k\ < 1. 


(6.18) 


This shows (6.17) for k £ Yf~ with |fc| < 1. The proof for k £ Yff with |fc| < 1 is similar. 


For k £ I3, we have | e ±t$ (C’ fe )| = 1 and 


—R 3 {x,t 1 k)Ri{x,t,k) — s 2 —{R±(x,t,k) + s)e 
{R 3 {x, t, k) — s)e t$ ^ ,fe ) 0 


u(x,t , k ) 
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Hence, by (6.10) 


\(u(x,t,k)) n| < \R 3 - s11i ? 4 1 + \s\\R.4 + s| < C\k\, 

\(u(x,t,k)) 2 i\ < |-R 3 ~ s\ < L\k\, \(u(x,t, k)) 12 \ < |A 4 + s\ < L\k\, 


for k G Y 3 and ( x,t ) G S. This shows (6.17) also for k G Y 3 and completes the proof of the 


claim. 


V 


t —» oo, (x, t) G S, 
t —> oo, ( x, t) G S, 


(6.19) 

( 6 . 20 ) 

( 6 . 21 ) 


Claim 2. We have 

INM,-)|| L oo(f) = 0(t~T>), 

\\w{ x ^r)\\ L 2 { f) = 

and 

\\w(x,t, -)IIli(t) = 0(t~%), oo, (x,t) G S, 

where the error terms are uniform with respect to x. 

Proof of Claim\f^ Let p = 2 or p = oo. Then 

l|w(M,-)||j> ( f) < II w(x,t, -)||Lp(r\y) + Wmoix^,-)- 1 - I\\ L p(\k\=i) + \\w{x,t, -)IUp(y)- (6-22) 

On T \ Y, the matrix ui is given by either w or niowm q 1 . Hence \\w(x,t, OIIl^^tvy) = 0(t~ s) 
and 

HwOMiOlU^rYY) < \j\\w(x,t, •)IU~(r\y)||w(a:,t, OlU^rvr) = 0(t~s) 
by the assumption ( |6Jl ) and the boundedness of nig. Moreover, by (|5.5[), 


m 0 (x,t,k) -/||iP(|fc|=i) = 


—x 


m [s,j—^_,(3t) 3 k,(3t) 3 k 0 ) -I 


-1 


LP(\k\=l) 


= o(t-s), 


(6.23) 


uniformly in x. By (6.1), we can choose a constant K > 0 such that 2kg < Kt. 1 / 3 for all 
(x,t) G S. Then, for |fc| > Kt _1 / 3 , we have 

C 

At\k\<- 

because equation ( 6. 14| ) implies 

C\k\ 


\w(x,t,k)\< ^—<Ct 3 , \k\> Kt 3 , k G Y, 


(6.24) 


\w(x,t,k)\ < C|fc|e - 4t|fc - fe ° |3 < C\k ^ k . fco|3e . 4t '„ fc < 


4f|fc — kg I 3 


4t|fc — kg\ 3 


c 


<^^<Ct—\ \k\>Ktk€Y+UY+, 

and a similar argument applies for k G Yf~ U Yf. On the other hand, if \k\ < AT -1 / 3 , equation 
(6.14) implies 


\w(x,t, k)\ < C\k\ < Ct 3 , \k\ < Kt 3 , k£Y. 


(6.25) 


Using equations (6.24) and (6.25), we can estimate the third term on the right-hand side of 


( 6 . 22 ) as follows 

IKz,t,-)|U~(y) = o(t~i), 


\w(x,t,-)\\ L 2( Y) < cj f (tTIu) +/ i l fc l 2 d l fe l 

y Jrn{|fc|>ift _ 3} ^*|«:| ' “'Yn{|fe|<ift _ 3} 


-| /*00 nt 3 

< ^ J i u~ A du + J u 2 du = 0(t~i). 


This proves (6.19) and (6.20). 
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In order to prove (6.21), we note that (6.14) implies 


\\w(x,t,-)\\ L i(x)<cf \w(x,t,k)\\dk\ + C f \w(x,t,k)\\dk\ 
JYiUY 2 Jy 3 

r°° 3 r k ° 

< C / (u + ko)e~ 4tu du + C / \u\du 

J 0 J -fcn 


0 J -k 0 

r (§) 2 , k 0 T(l)i 


~ C 6V% t t+C 22/3 


t 3 + Cfcg < Ct 3 , (i,t) eS. (6.26) 


This proves (6.21). 


V 


Let C denote the Cauchy operator associated with f: 

(OX*) = ^j/ S*. 


zec\r. 


We denote the nontangential boundary values of Cf from the left and right sides of f by C + f 
and C-f respectively. We dehne C-a, : L 2 (t) + L°°(t) — > L 2 (t) by C$f = C-(fw). 

Claim 3. There exists aT > 0 such that I — Cw( x ,t,-) € B(L 2 (Y)) is invertible for all (x,t) £ S 
with t > T. 

Proof of C/azm[^ The assumption that 5 r is uniformly bounded on L 2 (T) together with the 
Sokhotski-Plemelj formula C_ = |(—/ + Sf) show that sup( x ^g§ < oo. Thus, by 

(p9l), 


ll^v’|| i g(^ 2 (p)) Is C f ||d)||^ 00 (p) 0(t 3 ), t —> oo, (x,t) £ S, 

uniformly with respect to x. This proves the claim. 


(6.27) 

V 


In view of Claim |3j we may define the 2 x 2-matrix valued function ft(x, t, k) whenever t > T 


by 

// = / + (/ — Cu,) 1 Cwi s i + l 2 (t). 

Claim 4. The function fj,(x, t, k ) satisfies 

\\p( x ,t, •) — /|| I/ 2 ( f) = 0(t~i), f-t oo, (x,t) £ S, 
where the error term is uniform with respect to x. 

Proof of Claim [|J We have 

" 1 


(6.28) 

(6.29) 


ll( J C w) lle(L3(f)) ^ 2^ ll C *llB (L 2 ( f)) - 1 ||A II 

j= 0 1 - ll ( -''*|| B (L 2 (f)) 


whenever ||Cu>||^(x, 2 (Tj) < 1- Using the bound sup^ eI ||C_|| B , i2 /f)) < oo and equation (6.27), we 
find 

IIA ~ ^1li 2 (f) = II (I ~ C-w) C"wl\\ L 2 (f) — II {I ~ &w) IIg(L 2 (f)) 11^— 0^)11L 2 (f) 

<- . (r) — < C\\w\\ L2(t , 

1 — l|£ii|lB(L 2 (f)) 


for all (x,t) £ S with t large enough. In view of (6.20), this gives (6.29). 


V 


Claim 5. There exists a unique solution rh £ I + E 2 (C \ f) of the L 2 -RH problem (6.13) 
whenever t > T. This solution is given by 


1 f ds 

m(x,t,k) = I+ C(p,w) = I + -—: / p.(x,t, s)w(x,t, s) 


2iri 


- k' 


(6.30) 
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Proof of Claim, [5| Uniqueness follows since detD = 1. Moreover, equation ( |6.28[ ) implies that 
ft — I = Cwfh- Hence, by the theory of L 2 -RH problems, rh = I + C(fiw) satisfies the L 2 -RH 
problem (16.131). V 


Claim 6. For each point ( x,t ) £ S with t > T, the nontangential limit of k(m(x,t 1 k) — I) as 
k —> oo exists and is given by 

^ If 

lim k(m(x,t,k) — I) = — -— : / jj,(x,t,k)w(x,t,k)dk. 
k—>oo 2m Jf< 

Proof of C/aim[c| Fix (x,t) £ S and let W a ^ = {a < argfc < 6} be a nontangential sector 
at oo with respect to T. This means that there exists a S > 0 such that W a -s,b+8 does not 
intersect Tn{|z| > i?} whenever R > 0 is large enough. It follows that there exists a c > 0 such 
that |s — k\ > c(|s| + |fc|) for all s £ T and all k £ W a ,b with |fc| sufficiently large. Thus, since 
fiw £ T 1 (T), dominated convergence implies 

\s\\ds\ 

pb'uu)yx, o, s)us n urn ; \y/j,w)\x, i, S) \ ■ 

/r s ~ k J r 

\s\\ds\ 


z 

lim 

k—yoo 


kds f z r 

( p,w)(x,t,s )--+ / (fiw)(x,t,s)ds < lim / \(p,w)(x,t,s)\-. ,. 

S — k J f k-¥oo J p |s — k\ 

t Z f i/- a\/ M Islldsl 
< lim / \(pw)(x,t,s)\ 

fc->oo J p c(|s| + |fc|) 


= 0. 


The claim now follows from (6.30). 

Claim [6] implies 

z z 

lim k(m(x, t 1 k) — I) = lim k(m(x, t, k) — I) 


V 


k—yoo 


k—> oo 


= —-—: [ fi,(x,t,k)(mo(x,t,k) 1 — I)dk — —— : / fi(x,t,k)w(x 1 t,k)dk. (6.31) 

2tt* V|fc|=i 2 tt* J t 


By (5.5) 


mo(x, t,k) 1 = m z ( s, ——, (3 t) 3 k, (3 1) 3 ko 


(3 t)i 

= / _B(M) +0 


-l 


where B(x,t) is defined by 


B(x, t) 


\ 

2 lu' 


(3t)zk \(3t) 3 k 2 J ’ 

V(3t)3 / 


t 3 \k\ —> oo, 


(6.32) 


(f=fr;s,0,-a) 

V (3 1) 3 / 


Using (6.23), ( 6.29[ ), and ( |6.32 ), we find 
/ p(x,t, k)(mo(x,t, k)^ 1 — I)dk 

J\k\=l 

= / (mo(x, t, k)^ 1 — I)dk + / (/t(x, t, k) — I)(mo(x, t, /c) _1 — I)dk 

J\k\ — 1 J\k\=l 


i\k\=: 
B(x, t) 


J k i ']7+ 0 ( t + °(ll£(Mr)-Al£2 ( f)||(m 0 (M,') 1 - -0IU* ( |fc|=i)) 


(3 t)l J\k\=i k 

= _ + o(t-f), t->-oc, (i,f)e S, 

(3t) 3 

uniformly with respect to x. On the other hand, 
fi((, t, k)w((, t, k)dk 


(6.33) 


< IIA - AU 2 (r)IH| L 2 (r) + |H| l i ( t). 
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By (6.20) and (6.29), both ||io||L 2 (r) a nd ||A — -^l|L 2 (r) a re Q( t 3 ) as t —> 00 . Moreover, the L 1 
norm of w is 0(t~ 3 ) on f \ Y by (6.71. Since equations (6.12) and (5.6) imply that mo + (x,t,k) 


and nio-(x,t,k) are uniformly bounded for ( x,t ) £ S and k £ Tfl {|fc| < 1}, it follows that the 
L l norm of w also is 0(t~ 3) on T \ Y. On the other hand, the L 1 norm of w is 0(t _ s) on Y 
by (6.21). Thus ||w||ii(r) = 0(t~i). We infer that 


fi(x, t , l)w(x, t, l)dl 


= 0(t 3 ), t —> 00 , (x,t) £ S, 


uniformly with respect to x. Equations (6.31|), (|6.33[), and (6.34) give (6.11). 


(6.34) 

□ 


7. Asymptotics in the self-similar sector 
The purpose of this section is to prove Theorem |7.1| which describes the asymptotic behavior 


of quarter plane solutions of the rnKdV equation in the self-similar sector (|1.4|). The proof relies 
on an application of the nonlinear steepest descent result of Theorem 6.1 


Theorem 7.1 (Asymptotics in the self-similar sector). Let r £ C n (M) and h £ C 5 (D 2 ) be two 


functions which satisfy the assumptions of Theorem 4-1 Given N > 0, let § = {t > 1, 0 < x < 


AT 1 / 3 } be the self-similar sector defined in (1-4)- Let s = ir (Ob 

Then there exists a T > 0 such that the L Z -RH problem (2.9) has a unique solution and the 
limit in (2.1C\ ) exists for all ( x , t) £ S with t > T. Moreover, the function u(x, t) defined by 
(2.10) satisfies 


i(x,t) = 


' _ x 

.(3i) 1/3 ’ 


,0,-s) 


(3i)a 


+ 0(t 3 ), t 00 , (x,t) £ S, 


(7.1) 


where the error term is uniform with respect to x, and u p (■’, s , 0, —s) denotes the solution of the 
Painleve II equation (5.4) corresponding to (si,S 2 ,S 3 ) = (s,0,— s) according to (A. 6). 


Before presenting the proof of Theorem |7.1[ we give a corollary. 


Theorem 7.1 applies whenever the functions r(k) and h(k) satisfy the stated assumptions, 
regardless of whether these functions derive from a quarter plane solution or not. If we do 
assume that r(k) and h(k) are defined in terms of some initial and boundary values uq(x) and 
{gj (t )}q via (2.3), we obtain the following corollary. 


Corollary 7.2 (Asymptotics for quarter plane solutions). Suppose uq, go, gi, 32 satisfy (2-4) 
with n = 11 and m = 4, that the associated spectral functions satisfy the global relation (2.5), 
and that the homogeneous L 2 -RH problem associated with (2.9) has only the trivial solution for 
each (x,t) £ [0, 00 ) x [0,oo). Then the function u{x,t) defined by (2.10) is a classical solution of 


1.1) in the quarter plane with initial data u(x, 0) = Uq(x) and boundary values d 3 x u{ 0, t ) = gj(t), 
j = 0,1,2. Moreover, in the self-similar sector, the solution u(x,t) satisfies 


i(x,t) = 0(t 3 ), t —> 00 , (x,t) £ 


uniformly with respect to x. 


Proof of Corollary \7.2\ Combining Theorem |2.1| and Theorem |7.1[ we infer the existence of a 
classical quarter plane solution u{x,t) which satisfies the asymptotic formula (7.1) in the self¬ 
similar sector. The formula <Z3> being uniformly valid as x —> 0, can only be consistent with 
the assumed decay of u{ 0, t) provided that the coefficient u p vanishes at x = 0. In fact, it turns 
out that the global relation forces the parameter s to vanish, which in turn implies that u p 
vanishes identically. To see this, we first note that the symmetries 


a(k) = a(—k), b(k) = b(—k), A(k) = A(—k), B(k ) = B{—k), 


imply that the four values a(0), 6(0), A(0), -B(O) are real. The global relation (2.51 and the unit 


determinant relation a(k)a{k ) — b(k)b(k) = 1 then show that 6,(0) = —6(0)/o(0), which gives 
s = ir( 0) = 0. For s = 0, the jump matrices in the Painleve II RH problem ( A.3|) reduce to 
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the identity matrix, so the solution u p (y ; 0,0, 0) corresponding to s = 0 is identically zero. This 
completes the proof of the corollary. □ 


Proof of Theorem 7.1 Just like the proof of Theorem |4.f| the proof consists of seven steps. 


Step 1: Introduce an analytic approximation of h(k). Lemma 4.4 applies also for 
(x, t ) £ S. Thus we find a decomposition 

h(k) = h a (t,k) + h r (t,k), k £ dD 1 , 
where h a and h r have the properties listed in Lemma |4.4| 

Step 2: Deform. Proceeding as in the proof of Theorem |4.1 


we note that M satisfies the 
L 2 -RH problem (2.9) iff the function M h (x,t,k ) defined by (4.12) satisfies the L 2 -RH problem 


( 4.13D . 

Step 3: Conjugate. As in the proof of Theorem |4.1[ we conjugate the RH problem to 
arrive at an appropriate factorization. However, instead of introducing <5(cj’, k ) by (4.15), we now 
define a function S(k) by the equation obtained from (4.151 by replacing fco with zero (recall 
that 4 -> 0 in the self-similar sector). Thus, let 

^( fc)- 1 0 

0 S(k) y 


A(k) = 


where 

5(k) = e -^ ln(1 - |r(s)|2) A, k £ <C\R. 
The function 5 satisfies the following jump condition across the real axis: 

M *) 


{ +W = TA| rWr 


k £ 


(7.2) 

(7.3) 


The symmetry r{k) = r(—k) implies that A obeys the symmetries (4.18). 
Lemma 7.3. The 2 x 2-matrix valued function A (k) satisfies 

A, A" 1 £ I + E 2 (C \ R) n E °°(C \ R). 


Proof The proof is similar to, but easier than, the proof of Lemma 4.5 


V 


By Lemma |7.3[ M h (x,t, k) satisfies the L 2 -RH problem (4.131 iff the function M defined by 
M(x,t,k) = e~^ a3 M h (x,t,k)A(k)e^ a3 
satisfies the L 2 -RH problem 

|M(z,V) £ I + E 2 (C \ S), 
yM + (x,t, k) = M-(x,t,k)J(x,t,k) for a.e. k £ S, 


(7.4) 

(7.5) 


where 


J(x,t,k ) = e " a3 A_ 1 (k)J h (x,t,k)A + (k)e" a3 

^iS{k)~ 2 h r (t,k)e t ^^’ k ' ) 1^ ' 11 

= < f !-k(^)| 2 i8-(k) 2 l J\r k { l ) \2 e ~ mck) ^ 

Mo- 1 ® (C,fc)\ 

k £ 9D 4 , 


k £ 


(7.6) 


1 —i5(k) 2 h r (t,k)e 

,0 1 


and we recall that E = dD\ U dD± U R is the contour shown in Figure [Tj The upshot of the 
above conjugation is that we can now factorize the jump across R as follows: 


J = 1 B[. 


k £ R. 


(7.7) 
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Figure 11. The contour T = £ U Y\ U Y 2 and the open sets V\ and V 2 . 


where 

B 1 (i<L|_(fc) _2 r 1 (fc)e t<s ^’ fe ) 1 

and the functions fi(fc) and r 4 (fc) are defined by 

r(k) 


B„ = 


1 —i5-{k) 2 r/i{k)e 

0 1 


n(k) = 


1 — r(k)r(k) 


r 4 (/c) = 


r(k) 


1 — r(k)r(k) 


The factors e ± 4 a3 in (7.41 are inserted for later convenience (more precisely, to ensure that 
(6.10) holds with s pure imaginary). 


Step 4: Introduce analytic approximations of r\(k) and r 4 (/c). Lemma 4.6 applies 
also for (x, t ) G S. Thus we find decompositions 

rj(k) = rj, a (x,t,k) + rj tr (x,t,k), j = 1,4, k G R, 


where r 7jQ and r+ r have the properties listed in Lemma 4.6 


Step 5: Deform again. Let T = EU YiU Y 2 and, for j = 1, 2, let Vj denote the open region 
between R and Y), see Figure 11 Write 

Bi — Bi r Bi a , B u — B u ^ r B uai 

where {B^ a , B U}a } and {7?/ jr , B U}r } denote the matrices {Bi,B u } with (n(fc), r 4 (fc)} replaced 
with {r+ a (fc), r 4>a (fc)} and |ri jr (fc),r 4iT .(fc)}, respectively. The estimate (4.29) implies that 

IB lta (x,t, •), B^(x,t,-)ei+ (E 2 n £°°)(Vi), 

\B u , a {x, t, •), B-^x, t, ■) 6 / + (E 2 n E°°)(V 2 ), 


for each (x, t) G S. Hence M satisfies the L 2 -RH problem (7.51 iff the function m(x, t , k) defined 
by 


M(x,t, k)B[ a (x, t, fc) -1 , 
m(x,t,k) = <( M{x,t,k)B u ^ a (x,t,k)~ 1 , 
M(x, t, k), 


k £\ 1 , 
k G V 2 , 
elsewhere, 


satisfies the L 2 -RH problem 


m(x, t, •) G I + E 2 (C \ T), 

m+(x, t, fc) = m_(x, t, k)v(x, t, k ) for a.e. k G T, 


(7.8) 
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Figure 12. 


The contours T and F satisfy the properties (n)-(r4) of Theorem 


6.1 


where, in view of (4.24) and (4.25), the jump matrix v is given by 

( 1 


v = 


0 


Ru,a 


yiS(k) 2 r ly a(x,t, k)e mc ” k) 1 )' 
(1 — i5(k) 2 r^ a (x, t, 


1 °'l 

iS(k) 2 h r (t,k)e t ^^ ,k ^ 1 J 

1 —i5(k) 2 h r (t,k)e~ t '^^’ k ^ 

.0 1 


k £ Y\, 

k £ Y 2 . 

k £ Y 3 , 
k £ R\Y 3 , 

k £ dDi, 
k £ dD\. 


(7.9) 


Step 6: Apply Theorem |6.1[ We verify that Theorem |6.1| can be applied to the contour 
r and the jump matrix v with s = i(L|_(0) -2 ri(0) and 

Ri(x, t, k) = i5{k)~ 2 ri^ a (x, t , k), R 2 {x, t, k) = — i5{k) 2 r^ a {x 1 1, k), 

R 3 (x,t,k) = iS + (k)~ 2 ri(k), R±(x,t,k) = —z<5_(fc) 2 r 4 (fc). 


The contours T and F are shown in Figure 12 The conditions (IT)-(r4) are clearly satisfied. 
Since the contour fcp T is independent of (x,t), a scaling argument shows that ||»SrII b(^ 2 (T)) ' s 
independent of (x,t) £ S. Moreover, the Cauchy singular operator 5{|fc| = i} associated with the 
unit circle is bounded. Hence a decomposition argument shows that 5p is uniformly bounded on 
L 2 (F) (cf. the proof of Theorem 4.23 in [9]). Equation (6.4) follows from (7.9) and the estimates 
in Lemmas 4.4 and 4.6 Clearly detz; = 1. The symmetry condition (6.6) follows from (7.9) and 


the symmetries of h r , rj ta , rj, r , and 5. 

We next verify (6.7). By part (d) of Lemma 


4.6 


the L 1 and L°° norms of w are Oft 2 ) on 
the L 1 and L°° norms of 


4.4 


1 \ Y 3 uniformly with respect to (x, t) £ S. By part (c) of Lemma 
1 are 0(t~ a) on dD 3 U dD± uniformly with respect to (x,t) £ S. This verifies (6.7). 


It remains to verify (6.10). Since r(fc) = r(—k), we have r(0) £ R. Moreover, the symmetries 
6(k) = 6(—k) = S(— fc) _1 imply <5_|_(0) = <5 + (0) = 5_(0) _1 . In particular, s £ zR. It also follows 
that 


Rj(x,t, 0) = 


*M0) 


-2 


r(0) 


i-k(o )| 2 


= s, 


('QI2 r(0) _ _ _ 

l0 -\o) 1 —|r(0)| 2 


3 = 3, 
3 = 4 . 
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To see that |s| < 1, we write 5(k) = e x ^ k \ where x(^) = / R ln(l — l 7 ’( s )| 2 )- The 
Sokhotski-Plemelj theorem implies 


x+(°) = M 1 ~ l r ( s )l 2 )"7 - ^M 1 - IK°)I 2 ) 


ds 


(7.10) 


Hence 


|<5+(0)|- 2 = g- 2Re x+(°) = 1 _ | r ( 0 )| 2 , 

which implies |s| = |r(0)|. The assumption sup fegR \r(k)\ < 1 now gives |s| < 1. 

Finally, we prove ( 6. 10| ) for Ri and R 3 ; the proofs for R 2 and R 4 are similar. Using the 
expression (7. 10[) for x+(0), we obtain, for Imfc > 0, 


X(k) -X+(0) =- - |r(s)| 2 ) - l^j ds + M 1 “ l r (0)| 2 ) 




1 


In 


/l-|r( S )| 2 \ 

ds 

'1 k 

f ds 

Vi-k(0)| 2 J 

s~k 

2 27ri J 

r s{s - k)_ 


ln(l ~ l r (0)| 2 )- 


The expression within square brackets vanishes. Hence 

'1 - |r(s)| 2 


X(k) - X+(0) = ~J~ I ~ ln ( 7 

2ni J R s VI 


— |r(0)| 2 J s — k 1 


ds 


ImA’ > 0. 


The function ln( i_J^(ojj 2 ) belongs to i7 1 (K); thus (cf. (4.21)) 


\x{k) — X+(0)| < C\k\, Im k > 0. 


The inequality (4.22) now yields 

|1 _ e 2(x(fc)-x+(0))| < 2\x(k) - x + me 2Mk) ~ x+m < C\ X {k) - x+(0)| < C|fc| (7.11) 


k G Y+, (7.12) 


whenever InrA; > 0. On the other hand, by (4.28), (4.29), and Lemma 7.3 
\Ri(x,t, k) - s\ < |<5(A:)| _2 |ri ja (a:, A, k) - ri(fe 0 )| + |<5(fc)| —2 1(fc 0 ) - ri(0)| 

+ ~ ^+(0) _2 ||o(0)| 

< A: 0 |e4 |Re<I ' (c - fe)l +C|A:o| +C|l-e 2(x(fc) -x+(0))| ; 

where we recall that Y+ and Y ) - den ote t he parts of Yi that lie in the right and left half-planes, 
respectively. Equations (7.11) and (7.12) prove (6.10) for Ri(x,t,k) with k G Y +; the proof 
when k G Y{~ is analogous. The estimate (7.11) also implies 

I Ra(x,t,k) -s\< \ 6 + (k)\- 2 \ ri (k) - n(0)| + IMfc)- 2 - 5 + (0)- 2 ||n(0)| 

< C\k\ +C\l- e 2 (x+(O-x+(0))| < k e y 3i 


which proves (6.10) also for R 3 . 


Step 7: Find asymptotics. As a consequence of Theorem 6.1 the L 2 -RH problem (7.8) 


for m has a unique solution and the limit in (6.11) exists for all (x,t) G S with t sufficiently 
large. Since the RH problems for M and m are equivalent, the L 2 -RH problem (2.9) for M also 
has a unique solution for all sufficiently large t. Moreover, the limit (2.10) defining u(x,t) exists 
whenever t is large enough and 

z z 

u(x,t) = —2 i lim (kM(x,t,k))i 2 = 2 lim {km{x, A, k))\ 2 - 


k—>oo 


k—> 00 


Equation ( 6.11| ) of Theorem |6.1| then yields 

' P d3t)i/ 3;s >°’ s ) 


u{x, t) = 


(3 A) 3 


+ 0(t 3 ), t —> 00 , (i,() G S, 


where the error term is uniform with respect to x. This proves (7.1) and completes the proof of 
Theorem 17.11 □ 
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8. Concluding remarks 


Good asymptotic formulas provide one of the most efficient ways to obtain both quantita¬ 
tive and qualitative information on the solution of a PDE. Since the solution of an integrable 
equation can be expressed in terms of the solution of a matrix Riemann-Hilbert problem, the 
nonlinear steepest descent method introduced by Deift and Zhou m provides an excellent 
tool for establishing such formulas. For initial value problems, this method has been utilized 
extensively since its introduction in 1993. 

In this paper, we have considered the derivation of asymptotic formulas for initial- boundary 
value problems via the nonlinear steepest descent method. In particular, we have established 
asymptotic formulas valid in the similarity and self-similar sectors for the mKdV equation in 
the quarter plane. A special emphasis has been placed on giving complete proofs and on the 
derivation of rigorous and uniform error estimates. 

Previous studies addressing the question of asymptotics for initial-boundary value problems 
for integrable equations include oi nun [g. I n pa rticular, it was shown in 0 that the 
quarter plane solution u{x, t) of the mKdV equation (1.11 satisfies u(x, t) = 0(t -1 / 2 ) as (x, t) —> 


oo in a subsector of the similarity sector. Our results give the explicit form of the leading order 
term in this formula. We have also demonstrated that the boundary values have a large effect on 
the asymptotic behavior of the solution in the self-similar region. In particular, if the boundary 
values have sufficient decay as t —> oo, then we have showed that the solution is uniformly 
of 0(t~ 2/3 ) throughout the self-similar sector. In contrast, in the absence of a boundary, the 
solution has a leading order term of order 0(t -1 / 3 ) in this sector. 

Appendix A. The RH problem associated with Painleve II 


We give a number of results related to the Painleve II equation (5.4). We refer to [T3] for full 
proofs. Equation (5.4) is the compatibility condition 

Ay -U x + [A, U]= 0 

of the Lax pair 

(tf* = A*, 


'Ey = C/T, 


(A.l) 


where 


A(y, z ) = — *(4,2 2 + y + 2u 2 )a 3 — 4zua 2 — 2u y a\, 
U(y, z) = -iza 3 - ua 2 , 


and {cTj} 3 denote the standard Pauli matrices. The linear system (A.l) has an irregular singular- 
point at z = oo. The theory of ODEs with singular points implies that there exist solutions 
{T„}f of (A.l) such that 

*n(y,z) = (i + 0(z- 1 ))e- i ( yz +^ 3 , *-> 00 , zea n , 

where the sectors I ~l n are given by 


(A.2) 


— \ z g c 


ar g~ € ( |(n - 2), 


n = 1,... 6. 


The asymptotics in (A.2) is uniform as z —> oo within any closed subwedge of Q n . The T„’s are 


entire functions of z which are related by 

= 4'nW 

where the constant Stokes matrices {<S n }f have the triangular form 


^ n+ i = ^„5„, n = 1,..., 6, 


S n = 


1 0 

Sn 1 


n odd; 


S n = 


1 S r 
0 1 


n even. 


The symmetries 


A(y,z) = -<j 2 A(y,-z)a 2 , U(y,z) = cr 2 U(y,-z)a 2 
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Figure 13. The contour T = U® =i rj. 


imply that 


Sn +3 — i- 6 . S n _|_3 — 5 n , H — 1,2,3. 


Let T = U® = 1 T, where the rays {Tj}f are defined by (see Figure 13) 


r 3 = Uec 


7T 7T (n — 1) 

aT ^ ~ = (£ 3 


n = 1 ,..., 6 . 


Assume that these rays are oriented away from the origin. Then the sectionally analytic function 
m p (y,z) defined by 

' 7r(n — 2) 7 rn' 


m p (y,z) = ^ n (y,z)e^ z+ i z argz £ 


6 


6 


satisfies the RH problem 

j rn p (y , z) = m p (y, z)e~ l( - yz+ ^ z3 ' >a3 S n , z€f 


m p {y,z) = I + 0(z 1 ), 


n, n = 1 ,... ,6 
oo. 


(A.3) 


(A.4) 


The cyclicity condition S 1 S 2 S 3 S 4 S 5 S 6 = I reduces to the condition 

si - s 2 + s 3 + S1S2S3 = 0. 

It can be shown (see Theorem 3.4 and Theorem 4.2 in [T3j) that given a set of complex 
constants S = {si, s 2 , S3} satisfying (A.4), there exists a countable set Yg = {yj = yj(S)}jT 1 C 
C with yj —> 00 as j —> 00 , such that the RH problem (A.3) has a unique classical solution^ 
m p (y, z) for each y £ C\Yg. The functions m p (y, z), where m p denotes the restriction of mr 
to argz £ ( 7r ( 2 "~ 3 ) ; 7r ^ 2 ”~ 1 ^ ), admit analytic continuations to (C \ I 5 ) x C. As z —> 00 , m p 
satisfies 

u p {y\si ,s 2 ,s 3 ) 


m (y,z) = I + — 


0 


2 z \u p (y,si,s 2 ,s 3 ) 


0 


+ 0{z ), z -)• 00 , y $ Y s , 

(A.5) 


^A classical solution of (A.3) is a 2 x 2-matrix valued function m p (y , •) such that (i) for each component D 
of C \ r, the restri ction of rrPzo D is analytic in D and has a continuous extension to D , (ii) m p satisfies the 
jump condition in (A.3), and (Hi) m p (y,z) = I + 0(z —1 ) uniformly as z —» oo, 2 GC. 
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where u p (•; si, s 2 , S 3 ) satisfies the Painleve II equation (5.4). The asymptotics in (A.5) is uniform 
in the following sense: For any compact subset K C C \ I 5 , there exists a C > 0 such that 


1 


(y’Z)-!- a - 


0 


M p (y;si,s 2 ,s 3 ) 


2 z \u p (y;si,s 2 ,s 3 ) 


< C\i 


\z\ >1, y £ K. 


Moreover, the map S 1 


’(•; Si, s 2 , s 3 ) is a bijection 


{(si, s 2 , S 3 ) G CP | si - s 2 + s 3 + sis 2 s 3 = 0} -)■ {solutions of (5.4)} 


(A.6) 


and the set Y$ is the set of poles of u p {-\ Si, s 2 , s 3 ), see Corollary 4.4 in [Till . 

We will be interested in solutions u of the Painleve II equation which are real-valued on 
To this end, we note that (p. 158 in |13| ) 


(y; si, s 2 , s 3 ) = u p {y; si, s 2 , s 3 ) 


if and only if the s n ’s satisfy 


S3 Si, 


S 2 = 52- 
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